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These notes provide a pedagogical introduction to the theory of non-leptonic heavy- 
meson decays recently proposed by Beneke, Buchalla, Sachrajda and myself. We 
provide a rigorous basis for factorization for a large class of non-leptonic two-body 
B-meson decays in the heavy-quark limit. The resulting factorization formula 
incorporates elements of the naive factorization approach and the hard-scattering 
approach, and allows us to compute systematically radiative ( "non-factorizable" ) 
corrections to naive factorization for decays such as B — > Diz and B — > irn. 



1 Introduction 

Non-leptonic two-body decays of B mesons, although simple as far as the 
underlying weak decay of the b quark is concerned, are complicated on account 
of strong-interaction effects. If these effects could be computed, this would 
enhance tremendously our ability to uncover the origin of CP violation in 
weak interactions from data on a variety of such decays being collected at the 
B factories. In these lecture, I review recent progress towards a systematic 
analysis of weak heavy-meson decays into two energctk^mesons based on the 
factorization properties of decay amplitudes in QCDtrEI. My discussion will 
follow very closely the detailed account of this approach given inB. (We have 
worked so hard on this paper that any attempt to improve on it were bound to 
fail and leave the author in despair.) Much of the credit for these notes belongs 
to my collaborators Martin Beneke, Gerhard Buchalla, and Chris Sachrajda. 

As in the classic analysis of semi-leptonic B — > D transitions our ar- 
guments make extensive use of the fact that the b quark is heavy compared to 
the intrinsic scale of strong interactions. This allows us to deduce that non- 
leptonic decay amplitudes in the heavy-quark limit have a simple structure. 
The arguments to reach this conclusion, however, are quite different from those 
used for semi-leptonic decays, since for non-leptonic decays a large momentum 
is transferred to at least one of the final-state mesons. The results of our work 
justify naive factorization of four fermion operators for many, but not all, non- 
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lcptonic decays and imply that corrections termed "non-factorizable" , which 
up to now have been thought to be intractable, can be calculated rigorously if 
the mass of the decaying quark is large enough. This leads to a large number 
of predictions for CP-violating B decays in the heavy-quark limit, for which 
measurements will soon become available. 

Weak decays of heavy mesons involve three fundamental scales, the weak- 
interaction scale Myy, the 6-quark mass mj, and the QCD scale Aqcd, which 
are strongly ordered: Mw ^> nib Aqcd- The underlying weak decay be- 
ing computable, all theoretical work concerns strong-interaction corrections. 
QCD effects involving virtualities above the scale m& are well understood. 
They renormalize the coefficients of local operators Oi in the effective weak 
HamiltonianEl, so that the amplitude for the decay B — > M\Mi is given by 

A(B -> MM = ^| A * CiO-O {M 1 M 2 \O i { l j)\B) , (1) 

where each term in the sum is the product of a Cabibbo-Kobayashi-Maskawa 
(CKM) factor A^, a coefficient function C,(/^), which incorporates strong- 
interaction effects above the scale fi ~ m?,, and a matrix element of an operator 
Oi. The difficult theoretical problem is to compute these matrix elements or, 
at least, to reduce them to simpler non-perturbative objects. 

A variety of treatments of this problem exist, which rely on assumptions 
of some sort. Here we identify two somewhat contrary lines of approach. The 
first one, which we shall call "naive factorization" , replaces the matrix element 
of a four-fermion operatoiiin a heavy-quark decay by the product of the matrix 
elements of two currents En3, e.g. 

{D+n-\(cb) v „ A {du) V - A \B d ) -> (7r-|(duy_A|0) (D + \(cb)v- A \B d ) . (2) 

This assumes that the exchange of "non-factorizable" gluons between the ir~ 
and the (Bd D + ) system can be neglected if the virtuality of the gluons is below 
fi ~ nib- The non-leptonic decay amplitude then reduces to the product of a 
form factor and a decay constant. This assumption is in general not justified, 
except in the limit of a large number of colours in some cases. It deprives the 
amplitude of any physical mechanism that could account for rescattering in the 
final state. "Non-factorizable" radiative corrections must also exist, because 
the scale dependence of the two sides of (|j) is different. Since such corrections 
at scales larger than [i are taken into account in deriving the effective weak 
Hamiltonian, it appears rather arbitrary to leave them out below the scale y,. 
Various generalizations of the naive factorization approach have been proposed, 
which include new parameters that account for non-factorizable corrections. 
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In their most general form, these generalizations have nothing to do with the 
original "factorization" ansatz, but amount to a general parameterization of the 
matrix elements. Such general parameterizations are exact, but at the price of 
introducing many unknown parameters and eliminating any theoretical input 
on strong-interaction dynamics. 

The second method used to study non-leptonic decays is the hard-scatter- 
ing approach, which assumes the dominance of hard gluon exchange. The decay 
amplitude is then expressed as a convolution of a hard-scattering factor with 
light-cone wave functions of the participating mesons, in analogy with more 
familiar applications of this method to hard exclusive reactions involving only 
light hadronsQ'EI. In many cases, the hard-scattering contribution represents 
the leading term in an expansion in powers of Aqcd /Qi where Q denotes the 
hard scale. However, the short-distance dominance of hard exclusive processes 
is not enforced kinematically and relies crucially on the properties of hadronic 
wave functions. There is an important difference between light mesons and 
heavy mesons in this regard, because the light quark in a heavy meson at rest 
naturally has a small momentum of order Aqcd, while for fast light mesons a 
configuration with a soft quark is suppressed by the endpoint behaviour of the 
meson wave function. As a consequence, the soft (or Feynman) mechanism is 
power suppressed for hard exclusive processes involving light mesons, but it is 
of leading power for heavy-meson decays. 

It is clear from this discussion that a satisfactory treatment should take 
into account soft contributions, but also allow us to compute corrections to 
naive factorization in a systematic way. It is not at all obvious that such a 
treatment would result in a predictive framework. We will show that this does 
indeed happen for most non-leptonic two-body B decays. Our main conclusion 
is that "non-factorizable" corrections are dominated by hard gluon exchange, 
while the soft effects that survive in the heavy-quark limit are confined to the 
(BMi) system, where Mi denotes the meson that picks up the spectator quark 
in the B meson. This result is expressed as a factorization formula, which is 
valid up to corrections suppressed by powers of AQcu/mb- At leading power, 
non-perturbative contributions are parameterized by the physical form factors 
for the B — > Mi transition and leading-twist light-cone distribution amplitudes 
of the mesons. Hard perturbative corrections can be computed systematically 
in a way similar to the hard-scattering approach. On the other hand, because 
the B — > Mi transition is parameterized by a form factor, we recover the result 
of naive factorization at lowest order in a s . 

An important implication of the factorization formula is that strong rescat- 
tering phases are either perturbative or power suppressed in Aqcd /m^. It is 
worth emphasizing that the decoupling of Mi occurs in the presence of soft 



3 



interactions in the (BMi) system. In other words, while strong-interaction ef- 
fects in the B — > M% transition are not confined to small transverse distances, 
the other meson M2 is predominantly produced as a compact object with small 
transverse extension. The decoupling of soft effects then follows from "colour 
transparency" . The colour-transparency argument for exclusive B decays has 
already been noted in the literatureEJO, but it has never been developed into 
a factorization formula that could be used to obtain quantitative predictions. 

The approach described inEra is general and applies to decays into a heavy 
and a light meson (such as B — > Dir) as well as to decays into two light mesons 
(such as B — > 7T7r). Factorization does not hold, however, for decays such as 
B — > ttD and B — > DD, in which the meson that does not pick up the spectator 
quark in the B meson is heavy. For the main part in these lectures, we will 
focus on the case of B — ► D^L decays (with L a light meson), for which 
the factorization formula takes its simplest form, and power counting will be 
relatively straightforward. Occasionally, we will point out what changes when 
we consider more complicated decays such as B — > irir. A detailed treatment 
of these processes can be found inE3. 

The outline of these notes is as follows: In Sect. || we state the factoriza- 
tion formula in its general form. In Sect. || we collect the physical arguments 
that lead to factorization and introduce our power-counting scheme. We show 
how light-cone distribution amplitudes enter, discuss the heavy-quark scaling 
of the B — > D form factor, and explain the cancellation of soft and collinear 
contributions in "non-factorizable" vertex corrections to non-leptonic decay 
amplitudes. We also comment on the implications of our results for final- 
state interactions in hadronic B decays. The cancellation of long-distance 
singularities is demonstrated in more detail in Sect. |], where we present the 
calculation of the hard-scattering functions at one-loop order for decays into a 
heavy and a light meson. Various sources of power-suppressed effects, which 
give corrections to the factorization formula, are discussed in Sect. ||. They in- 
clude hard-scattering contributions, weak annihilation, and contributions from 
multi-particle Fock states. We then point out some limitations of the factor- 
ization approach. In Sect. we consider the phenomenology of B — > D^L 
decays on the basis of the factorization formula and discuss various tests of 
our theoretical framework. We also examine to what extent a charm meson 
should be considered as heavy or light. Section || contains the conclusion. 

2 Statement of the factorization formula 

In this section we summarize the factorization formula for non-leptonic B de- 
cays. We introduce relevant terminology and definitions. 
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2.1 The idea of factorization 

In the context of non-leptonic decays the term "factorization" is usually ap- 
plied to the approximation of the matrix element of a four-fermion operator by 
the product of a form factor and a decay constant, as illustrated in (|^). Correc- 
tions to this approximation are called "non-factorizable" . We will refer to this 
approximation as "naive factorization" and use quotes on "non-factorizable" to 
avoid confusion with the (much less trivial) meaning of factorization in the con- 
text of hard processes in QCD. In the latter case, factorization refers to the sep- 
aration of long-distance contributions to the process from a short-distance part 
that depends only on the large scale nib. The short-distance part can be com- 
puted in an expansion in the strong coupling a s (rnb). The long-distance contri- 
butions must be computed non-perturbatively or determined experimentally. 
The advantage is that these non-perturbative parameters are often simpler 
in structure than the original quantity, or they are process independent. For 
example, factorization applied to hard processes in inclusive hadron-hadron 
collisions requires only parton distributions as non-perturbative inputs. Par- 
ton distributions are much simpler objects than the original matrix element 
with two hadrons in the initial state. On the other hand, factorization applied 
to the B — ► D form factor leads to a non-perturbative object (the "Isgur- 
Wise function"), which is still a function of the momentum transfer. However, 
the benefit here is that symmetries relate this function to other form factors. 
In the case of non-leptonic B decays, the simplification is primarily of the 
first kind (simpler structure). We call those effects non-factorizable (without 
quotes) which depend on the long-distance properties of the B meson and both 
final-state mesons combined. 

The factorization properties of non-leptonic decay amplitudes depend on 
the two-meson final state. We call a meson "light" if its mass m remains finite 
in the heavy-quark limit. A meson is called "heavy" if its mass scales with mb 
in the heavy-quark limit, such that m/m,b stays fixed. In principle, we could 
still have m Aqcd for a light meson. Charm mesons could be considered as 
light in this sense. However, unless otherwise mentioned, we assume that m is 
of order Aqcd for a light meson, and we consider charm mesons as heavy. In 
evaluating the scaling behaviour of the decay amplitudes, we assume that the 
energies of both final-state mesons (in the i3-meson rest frame) scale with m/, 
in the heavy-quark limit. 

2.2 The factorization formula 

We consider a generic weak decay B — > M\Mi in the heavy-quark limit and 
differentiate between decays into final states containing a heavy and a light 
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Figure 1: Graphical representation of the factorization formula. Only one of the two 
form-factor terms in (H) is shown for simplicity. 



meson or two light mesons. Our goal is to show that, up to power corrections 
of order AQCD/^bi the transition matrix element of an operator Oj in the 
effective weak Hamiltonian can be written as 

{M 1 M 2 \O i \B) = Y]if^ Ml (w|)/M 2 / duT^{u)® M M 

if M\ is heavy and M2 is light, 

{MiM 2 \Oi\B)=Y^Ff^{ml)f M2 f duT^u) * Ma (u) + {M 1 ^M 2) 

Jo 



+ /b/mJm 2 / dZdudvTf I (Z,u,v)&B(€)®M 1 (v)®M 2 ( u ) 
Jo 

if Mi and M 2 are both light. (3) 

Here F®^ M (m 2 ) denotes a B — ► M form factor evaluated at q 2 = m 2 , m 12 
are the light meson masses, and is the light-cone distribution amplitude 

for the quark-antiquark Fock state of the meson X. These non-perturbative 
quantities will be defined below. T^(u) and T* 1 (£,u,v) are hard-scattering 
functions, which are perturbatively calculable. The factorization formula in 
its general form is represented graphically in Fig. [j]. 

The second equation in (^) applies to decays into two light mesons, for 
which the spectator quark in the B meson (in the following simply referred 
to as the "spectator quark") can go to either of the final-state mesons. An 
example is the decay B~ — > n°K~. If the spectator quark can go only to 
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one of the final-state mesons, as for example in Bd — ► tt + K~, we call this 
meson Mi, and the second form- factor term on the right-hand side of (||) 
is absent. The formula simplifies when the spectator quark goes to a heavy 
meson (first equation in (||)), such as in Bd — > D + ir~. Then the second term 
in Fig. [l], which accounts for hard interactions with the spectator quark, can 
be dropped because it is power suppressed in the heavy-quark limit. In the 
opposite situation that the spectator quark goes to a light meson but the other 
meson is heavy, factorization does not hold, because the heavy meson is neither 
fast nor small and cannot be factorized from the B — ► M\ transition. Finally, 
notice that annihilation topologies do not appear in the factorization formula, 
since they do not contribute at leading order in the heavy-quark expansion. 

Any hard interaction costs a power of a s . As a consequence, the hard- 
spectator term in the second formula in (||) is absent at order a®. Since at this 
order the functions T^(u) are independent of u, the convolution integral results 
in the normalization of the meson distribution amplitude, and (^) reproduces 
naive factorization. The factorization formula allows us to compute radiative 
corrections to this result to all orders in a s ■ Further corrections are suppressed 
by powers of Aqcd/w!, in the heavy-quark limit. 

The significance and usefulness of the factorization formula stems from 
the fact that the non-perturbative quantities appearing on the right-hand side 
of the two equations in (||) are much simpler than the original non-leptonic 
matrix elements on the left-hand side. This is because they either reflect 
universal properties of a single meson (light-cone distribution amplitudes) or 
refer only to a B — > meson transition matrix element of a local current (form 
factors). While it is extremely difficult, if not impossible 113, to compute the 
original matrix element (M\M2\Oi\B) in lattice QCD, form factors and light- 
cone distribution amplitudes are already being computed in this way, although 
with significant systematic errors at present. Alternatively, form factors can 
be obtained using data on semi-lcptonic decays, and light-cone distribution 
amplitudes by comparison with other hard exclusive processes. 

After having presented the most general form of the factorization formula, 
we will from now on restrict ourselves to the case of heavy-light final states. 
Then the (simpler) first formula in (|^) applies, and only the first term shown 
in Fig. [I] is present at leading power. 

2.3 Definition of non-perturbative parameters 

The form factors F^^ M (q 2 ) in (^) arise in the decomposition of current matrix 
elements of the form (M(p')\qTb\B(p)), where T can be any irreducible Dirac 
matrix that appears after contraction of the hard subgraph to a local vertex 
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with respect to the B — > M transition. We will often refer to the matrix 
element of the vector current evaluated between a B meson and a pseudoscalar 
meson P, which is conventionally parameterized as 



(P(p')\q^b\B(p)) = F^(q 2 ) (p" +p'") 



" 2 <A ( 4 ) 
J g- 

where 5 = p — p', and .F^ - ^(0) = F^~* p (0) at zero momentum transfer. 
Note that we write (||) in terms of physical form factors. In principle, Fig. [I] 
could be looked upon in two different ways. We could suppose that the region 
represented by Fj accounts only for the soft contributions to the B — > Mi form 
factor. The hard contributions to the form factor would then be considered as 
part of TL (or as part of the second diagram) . Performing this split-up would 
require that one understands the factorization of hard and soft contributions to 
the form factor. If Mi is heavy, this amounts to matching the form factor onto 
a form factor defined in heavy-quark effective theorylij. However, for a light 
meson Mi the factorization of hard and soft contributions to the form factor 
is not yet completely understood. We bypass this problem by interpreting Fj 
as the physical form factor, including hard and soft contributions. This avoids 
the above problem, and in addition has the advantage that the physical form 
factors are directly related to measurable quantities. 

Light-cone distribution amplitudes play the same role for hard exclusive 
processes that parton distributions play for inclusive processes. As in the latter 
case, the leading-twist distribution amplitudes, which are the ones we need 
at leading power in the 1/mf, expansion, are given by two-particle operators 
with a certain hclicity structure. The hclicity structure is determined by the 
angular momentum of the meson and the fact that the spinor of an energetic 
quark has only two large components. The leading- twist light-cone distribution 
amplitudes for pseudoscalar mesons (P) and longitudinally polarized vector 
mesons (Vji) with flavour content (qq') are defined as 

(P(q)\q(y) a q'(x) \O) = ^ (rfys)^ J due l ^+^ $p{u,h), 
(V ll (q)\q(y) a q'(x)(l\0) = ~4p a [ due^+^^u^), (5) 



where (x — y) 2 = 0. We have suppressed the path-ordered exponentials that 
connect the two quark fields at different positions and make the light-cone 
operators gauge invariant. The equality sign is to be understood as "equal 
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up to higher-twist terms". It is also understood that the operators on the 
left-hand side are colour singlets. When convenient, we use the "bar" -notation 
U = 1 — u. The parameter fi is the renormalization scale of the light-cone 
operators on the left-hand side. The distribution amplitudes are normalized as 
J du$x(u, /i) = 1 with X = P, Vj| . One defines the asymptotic distribution 
amplitude as the limit in which the renormalization scale is sent to infinity. In 
this case 

$x(w,/i) P =°°6 M (l- W ). (6) 

The use of light-cone distribution amplitudes in non-leptonic B decays 
requires justification, which we will provide in Sects. || and |j. The decay 
amplitude for a B decay into a heavy-light final state is then calculated by 
assigning momenta uq and uq to the quark and antiquark in the outgoing light 
meson (with momentum q), writing down the on-shell amplitude in momentum 
space, and performing the replacement 

U aa (uq) T(u, . . -)al3,abVl3b(uq) — > — — / du$ P (u) (fa) /3a r(lt, . . .) a p,aa (7) 

4jv c J Q 

for pseudoscalars and, with obvious modifications, for vector mesons. (Even 
when working with light-cone distribution amplitudes it is not always justified 
to perform the collinear approximation on the external quark and antiquark 
lines right away. One may have to keep the transverse components of the 
quark and antiquark momenta until after some operations on the amplitude 
have been carried out. However, these subtleties do not concern calculations 
at leading-twist order.) 

3 Arguments for factorization 

In this section we provide the basic power-counting arguments that lead to 
the factorized structure shown in (^). We do so by analyzing qualitatively the 
hard, soft and collinear contributions to the simplest Feynman diagrams. 

3.1 Preliminaries and power counting 

For concreteness, we label the charm meson which picks up the spectator quark 
by M\ = D + and assign momentum p' to it. The light meson is labeled 
M% = 7r~ and assigned momentum q = En + , where E is the pion energy in 
the B rest frame, and n± = (1, 0, 0, ±1) are four- vectors on the light-cone. At 
leading power, we neglect the mass of the light meson. 

The simplest diagrams that we can draw for a non-leptonic decay ampli- 
tude assign a quark and antiquark to each meson. We choose the quark and 
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antiquark momenta in the pion as 



l q = uq + l ± + j-^ n_ , l q = uq-l± + n_ . (8) 

Note that q ^ l q + lq, but the off-shellness (i g + Z^) 2 is of the same order as the 
light meson mass, which we can neglect at leading power. A similar decompo- 
sition (with longitudinal momentum fraction v and transverse momentum l' ± ) 
is used for the charm meson. 

To prove the factorization formula (||) for the case of heavy-light final 
states, one has to show that: 

i) There is no leading (in powers of Aqcrj/rn^) contribution to the ampli- 
tude from the endpoint regions u ~ J^QGT>/ m b and u ~ AQcu/ m b- 

ii) One can set l± = in the amplitude (more generally, expand the ampli- 
tude in powers of l±) after collinear subtractions, which can be absorbed 
into the pion wave function. This, together with i), guarantees that the 
amplitude is legitimately expressed in terms of the light-cone distribution 
amplitudes of pion. 

hi) The leading contribution comes from v ~ (the region where 

the spectator quark enters the charm meson as a soft parton), which 
guarantees the absence of a hard spectator interaction term. 

iv) After subtraction of infrared contributions corresponding to the light- 
cone distribution amplitude and the form factor, the leading contribu- 
tions to the amplitude come only from internal lines with virtuality that 
scales with nib- 

v) Non-valence Fock states are non-leading. 

The requirement that after subtractions virtualities should be large is ob- 
vious to guarantee the infrared fmiteness of the hard-scattering functions Tfj. 
Let us comment on setting transverse momenta in the wave functions to zero 
and on endpoint contributions. Neglecting transverse momenta requires that 
we count them as order Aqcd when comparing terms of different magnitude 
in the scattering amplitude. This conforms to our intuition and the assump- 
tion of the parton model, that intrinsic transverse momenta are limited to 
hadronic scales. However, in QCD transverse momenta are not limited, but 
logarithmically distributed up to the hard scale. The important point is that 
contributions that violate the starting assumption of limited transverse mo- 
mentum can be absorbed into the universal light-cone distribution amplitudes. 
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The statement that transverse momenta can be counted of order Aqcd is to 
be understood after these subtractions have been performed. 

The second comment concerns endpoint contributions in the convolution 
integrals over longitudinal momentum fractions. These contributions are dan- 
gerous, because we may be able to demonstrate the infrared safety of the 
hard-scattering amplitude under assumption of generic u and independent of 
the shape of the meson distribution amplitude, but for « -> or u 1 a 
propagator that was assumed to be off-shell approaches the mass-shell. If such 
a contribution were of leading power, we would not expect the perturbative 
calculation of the hard-scattering functions to be reliable. 

Estimating endpoint contributions requires knowledge of the endpoint be- 
haviour of the light-cone distribution amplitude. Since it enters the factoriza- 
tion formula at a renormalization scale of order mj, we can use the asymptotic 
form (^) to estimate the endpoint contribution. (More generally, we only have 
to assume that the distribution amplitude at a given scale has the same end- 
point behaviour as the asymptotic amplitude. This is generally the case, unless 
there is a conspiracy of terms in the Gegenbauer expansion of the distribution 
amplitude. If such a conspiracy existed at some scale, it would be destroyed 
by evolving the distribution amplitude to a different scale.) We count a light- 
meson distribution amplitude as order Aqcd / mb i n the endpoint region (de- 
fined as the region the quark or antiquark momentum is of order Aqcd), and 
order 1 away from the endpoint, i.e. (for X = P, Vji) 



Note that the endpoint region has a size of order Aqcd/ n^b, so that the end- 
point suppression is ~ (AQco/wb) 2 . This suppression has to be weighted 
against potential enhancements of the partonic amplitude when one of the 
propagators approaches the mass shell. The counting for B mesons, or heavy 
mesons in general, is different. Naturally, the heavy quark carries almost all 
of the meson momentum, and hence we count 



The zero probability for a light spectator with momentum of order rrib must 
be understood as a boundary condition for the wave function renormalized at 
a scale much below m&. There is a small probability for hard fluctuations that 
transfer large momentum to the spectator. This "hard tail" is generated by 




1 ; generic u, 

Aqcd/w6; u, u ~ Aqcd/w6. 



(9) 




£ ~ AQ CD /m h , 



(10) 
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(a) (b) 

Figure 2: Leading contributions to the B — > D form factor in the hard-scattering approach. 
The dashed line represents the weak current. The two lines to the left belong to the B 
meson, the ones to the right to the recoiling charm meson. 

evolution of the wave function from a hadronic scale to a scale of order mj,. If 
we assume that the initial distribution at the hadronic scale falls sufficiently 
rapidly for £ 3> Aqcd/w?,, this remains true after evolution. We shall assume 
a sufficiently fast fall-off, so that, for the purposes of power counting, the 
probability that the spectator-quark momentum is of order nib can be set to 
zero. The same counting applies to the D meson. (Despite the fact that the 
charm meson has momentum of order mj, we do not need to distinguish the 
rest frames of B and D for the purpose of power counting, because the two 
frames are not connected by a parametrically large boost. In other words, the 
components of the spectator quark in the D meson are still of order Aqcd-) 

3.2 The B — > D form factor 

We now demonstrate that the B — ► D form factor receives a leading contribu- 
tion from soft gluon exchange. This implies that a non-leptonic decay cannot 
be treated completely in the hard-scattering picture, and so the form factor 
should enter the factorization formula as a non-perturbative quantity. 

Consider the diagrams shown in Fig. |^. When the exchanged gluon is 
hard the spectator quark in the final state has momentum of order mt,. But 
according to the counting rule ( |l0| ) this configuration has no overlap with the 
D-meson wave function. On the other hand, there is no suppression for soft 
gluons in Fig. ^. It follows that the dominant behaviour of the B — > D form 
factor in the heavy-quark limit is given by soft processes. 

Because of this argument, we can exploit the heavy-quark symmetries to 
determine how the form factor scales in the heavy-quark limit. The well-known 
result is that the form factor scales like a constant (modulo logarithms), since 
it is equal to one at zero velocity transfer and independent of as long 
as the Lorentz boost that connects the B and D rest frames is of order 1. 
The same conclusion follows from the power-counting rules for light-cone wave 
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Figure 3: Leading-order contribution to the hard-scattering kernels T?Au). The weak decay 
of the b quark through a four-fcrmion operator is represented by the black square. 

functions. To see this, we represent the form factor by an overlap integral of 
wave functions (not integrated over transverse momentum), 

F^ D (0)~ J *^* B frk x )* D (?®,k x ), (11) 

where £'(£) is fixed by kinematics, and we have set q 2 = for simplicity. The 
probability of finding the B meson in its valence Fock state is of order 1 in the 
heavy-quark limit, i.e. 

/^^IW£,^)| 2 ~i. (12) 

Counting k± ~ Aqcd and d£ ~ Aqcd /nib, we deduce that ^s^, fcjj ~ 
, we then obtain the scaling law F+~^ D (0) ~ 1, in 
agreement with the prediction of heavy-quark symmetry. 

The representation ( pi] ) of the form factor as an overlap of wave functions 
for the two-particle Fock state of the heavy meson is not rigorous, because 
there is no reason to assume that the contribution from higher Fock states 
with additional soft gluons is suppressed. The consistency with the estimate 
based on heavy-quark symmetry shows that these additional contributions are 
not larger than the two-particle contribution. 

3.3 Non-leptonic decay amplitudes 

We now turn to a qualitative discussion of the lowest-order and one-gluon 
exchange diagrams that could contribute to the hard-scattering kernels T^{u) 
in (||). In the figures which follow, the two lines directed upwards represent 
7r~, the lines on the left represent Bd, and the lines on the right represent D + . 

Lowest-order diagram 

There is a single diagram with no hard gluon interactions shown in Fig. ||. Ac- 
cording to (|To|) the spectator quark is soft, and since it does not undergo a hard 
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Figure 4: Diagrams at order a s that need not be calculated. 

interaction it is absorbed as a soft quark by the recoiling meson. This is evi- 
dently a contribution to the left-hand diagram of Fig. [|, involving the B — > D 
form factor. The hard subprocess in Fig. [| is just given by the insertion of 
a four-fermion operator, and hence it does not depend on the longitudinal 
momentum fraction u of the two quarks that form the emitted 7r~. Conse- 
quently, the lowest-order contribution to T^{u) in (j^) is independent of u, and 
the M-integral reduces to the normalization condition for the pion distribution 
amplitude. The result is, not surprisingly, that the factorization formula repro- 
duces the result of naive factorization if we neglect gluon exchange. Note that 
the physical picture underlying this lowest-order process is that the spectator 
quark (which is part of the B — ► D form factor) is soft. If this is the case, the 
hard-scattering approach misses the leading contribution to the non-leptonic 
decay amplitude. 

Putting together all factors relevant to power counting, we find that in 
the heavy-quark limit the decay amplitude for a decay into a heavy-light final 
state (in which the spectator quark is absorbed by the heavy meson) scales as 

A{B d -» D+7T-) ~ G F m 2 b F B ^ D {0) f n ~ G F m 2 b A QC d . (13) 

Other contributions must be compared with this scaling rule. 

Factorizable diagrams 

In order to justify naive factorization as the leading term in an expansion in a s 
and AQCD/^h, we must show that radiative corrections are either suppressed 
in one of these two parameters, or already contained in the definition of the 
form factor and the pion decay constant. Consider the graphs shown in Fig. |]. 
The first three diagrams are part of the form factor and do not contribute to 
the hard-scattering kernels. Since the first and third diagrams contain leading 
contributions from the region in which the gluon is soft, they should not be 
considered as corrections to Fig. [| However, this is of no consequence since 
these soft contributions are absorbed into the physical form factor. 

The fourth diagram in Fig. |] is also factorizable. In general, this graph 
would split into a hard contribution and a contribution to the evolution of the 
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Figure 5: "Non-factorizable" vertex corrections. 

pion distribution amplitude. However, as the leading-order diagram in Fig. |^ 
involves only the normalization integral of the pion distribution amplitude, the 
sum of the fourth diagram in Fig. |4| and the wave-function renormalization of 
the quarks in the emitted pion vanishes. In other words, these diagrams would 
renormalize the (ud) light-quark current, which however is conserved. 

"Non-factorizable" vertex corrections 

We now begin the analysis of "non-factorizable" diagrams, i.e. diagrams con- 
taining gluon exchanges that cannot be associated with the B->D form factor 
or the pion decay constant. At order a s , these diagrams can be divided into 
three groups: vertex corrections, hard spectator interactions, and annihilation 
diagrams. 

The vertex corrections shown in Fig. [5] viola±£L the naive factorization 
ansatz (||). One of the key observations made in Era is that these diagrams 
are calculable nonetheless. Let us summarize the argument here, postpon- 
ing the explicit evaluation of these diagrams to Sect. [|. The statement is 
that the vertex-correction diagrams form an order-a s contribution to the hard- 
scattering kernels T^(u). To demonstrate this, we have to show that: i) The 
transverse momentum of the quarks that form the pion can be neglected at 
leading power, i.e. the two momenta in (Q) can be approximated by uq and 
uq, respectively. This guarantees that only a convolution in the longitudinal 
momentum fraction u appears in the factorization formula, ii) The contribu- 
tion from the soft-gluon region and gluons collinear to the direction of the pion 
is power suppressed. In practice, this means that the sum of these diagrams 
cannot contain any infrared divergences at leading power in Aqcd/to6- 

Neither of the two conditions holds true for any of the four diagrams 
individually, as each of them separately contains collinear and infrared di- 
vergences. As will be shown in detail later, the infrared divergences cancel 
when one sums over the gluon attachments to the two quarks comprising the 
emission pion ((a+b), (c+d) in Fig. ||). This cancellation is a technical mani- 
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festation of Bjorken's colour-transparency argumenttL9, stating that soft gluon 
interactions with the emitted colour-singlet (ud) pair are suppressed because 
they interact with the colour dipole moment of the compact light-quark pair. 
Collinear divergences cancel after summing over gluon attachments to the b 
and c quark lines ((a+c), (b+d) in Fig. [5)). Thus the sum of the four dia- 
grams (a-d) involves only hard gluon exchange at leading power. Because the 
hard gluons transfer large momentum to the quarks that form the emission 
pion, the hard-scattering factor now results in a non-trivial convolution with 
the pion distribution amplitude. "Non-factorizable" contributions are there- 
fore non-universal, i.e. they depend on the quantum numbers of the final-state 
mesons. 

Note that the colour-transparency argument, and hence the cancellation of 
soft gluon effects, applies only if the (ud) pair is compact. This is not the case 
if the emitted pion is formed in a very asymmetric configuration, in which one 
of the quarks carries almost all of the pion momentum. Since the probability 
for forming a pion in such an endpoint configuration is of order (Aqcd /fni,) 2 , 
they could become important only if the hard-scattering amplitude favoured 
the production of these asymmetric pairs, i.e. if ~ 1/u 2 for u — > (or 
Tjj ~ 1/u 2 for u — * 1). However, we will see that such strong endpoint 
singularities in the hard-scattering amplitude do not occur. 

To complete the argument, we have to show that all other types of con- 
tributions to the non-leptonic decay amplitudes are power suppressed in the 
heavy-quark limit. This includes interactions with the spectator quark, weak 
annihilation graphs, and contributions from higher Fock components of the 
meson wave functions. This will be done in Sect. || In summary, then, for 
hadronic B decays into a light emitted and a heavy recoiling meson the first 
factorization formula in (|3|) holds. At order a s , the hard-scattering kernels 
Tlj(u) are computed from the diagrams shown in Figs. ^ and || Naive factor- 
ization follows when one neglects all corrections of order Aqcd / mh o,nd a s . 
The factorization formula allows us to compute systematically corrections to 
higher order in a s , but still neglects power corrections. 

3.4 Remarks on final-state interactions 

Some of the loop diagrams entering the calculation of the hard-scattering ker- 
nels have imaginary parts, which contribute to the strong rescattering phases. 
It follows from our discussion that these imaginary parts are of order a s or 
Aqcd/w;,. This demonstrates that strong phases vanish in the heavy-quark 
limit (unless the real parts of the amplitudes are also suppressed). Since this 
statement goes against the folklore that prevails from the present understand- 
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ing of this issue, and since the subject of final-state interactions (and of strong- 
interaction phases in particular) is of paramount importance for the interpre- 
tation of CP-violating observables, a few additional remarks are in order. 

Final-state interactions are usually discussed in terms of intermediate 
hadronic states. This is suggested by the unitarity relation (taking B — ► nir 
for definitcness) 



where n runs over all hadronic intermediate states. We can also interpret 
the sum in (JIJ) as extending over intermediate states of partons. The par- 
tonic interpretation is justified by the dominance of hard rescattering in the 
heavy-quark limit. In this limit, the number of physical intermediate states is 
arbitrarily large. We may then argue on the grounds of parton-hadron duality 
that their average is described well enough (up to Aqcd corrections, say) 
by a partonic calculation. This is the picture implied by (^]) . The hadronic lan- 
guage is in principle exact. However, the large number of intermediate states 
makes it intractable to observe systematic cancellations, which usually occur 
in an inclusive sum over hadronic intermediate states. 

A particular contribution to the right-hand side of ( |T4j ) is elastic rescatter- 
ing (n = 7T7r) . The energy dependence of the total elastic 7T7r-scattering cross 
section is governed by soft pomeron behaviour. Hence the strong-interaction 
phase of the B — > ttt: amplitude due to elastic rescattering alone increases 
slowly in the heavy-quark limitLJ. On general grounds, it is rather improbable 
that elastic rescattering gives an appropriate representation of the imaginary 
part of the decay amplitude in the heavy-quark limit. This expectation is 
also borne out in the framework of Regge behaviour, as discussed intJ, where 
the importance (in fact, dominance) of inelastic rescattering was emphasized. 
However, this discussion left open the possibility of soft rescattering phases that 
do not vanish in the heavy-quark limit, as well as the possibility of systematic 
cancellations, for which the Regge approach does not provide an appropriate 
theoretical framework. 

Eq. (||) implies that such systematic cancellations do occur in the sum over 
all intermediate states n. It is worth recalling that similar cancellations are 
not uncommon for hard processes. Consider the example of e + e _ -^hadrons 
at large energy q. While the production of any hadronic final state occurs 
on a time scale of order 1/Aqcd (and would lead to infrared divergences if 
we attempted to describe it using perturbation theory), the inclusive cross 
section given by the sum over all hadronic final states is described very well 
by a (qq) pair that lives over a short time scale of order 1/q. In close analogy, 
while each particular hadronic intermediate state n in (O) cannot be described 
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partonically, the sum over all intermediate states is accurately represented by 
a (qq) fluctuation of small transverse size of order Because the (qq) pair 

is small, the physical picture of rescattering is very different from elastic rnr 
scattering. 

In perturbation theory, the pomeron is associated with two-gluon ex- 
change. The analysis of two-loop contributions to the non-leptonic decay am- 
plitude inH shows that the soft and collinear cancellations that guarantee the 
partonic interpretation of rescattering extend to two-gluon exchange. Hence, 
the soft final-state interactions are again subleading as required by the validity 
of (||). As far as the hard rescattering contributions are concerned, two-gluon 
exchange plus ladder graphs between a compact (qq) pair with energy of or- 
der rrib and transverse size of order l/m& and the other pion does not lead 
to large logarithms, and hence there is no possibility to construct the (hard) 
pomeron. Note the difference with elastic vector-meson production through a 
virtual photon, which also involves a compact (qq) pair. However, in this case 
one considers s ^> Q 2 , where y/s is the photon-proton center-of-mass energy 
and Q the virtuality of the photon. This implies that the (qq) fluctuation is 
born long before it hits the proton. It is this difference of time scales, non- 
existent in non-leptonic B decays, that permits pomeron exchange in elastic 
vector-meson production in 7*p collisions. 



4 B — y Dir: Factorization at one-loop order 

We now present a more detailed treatment of the exclusive decays B c i — > 
DW+L-, where L is a light meson. We illustrate explicitly how factoriza- 
tion emerges at one-loop order and compute the hard-scattering kernels T^(u) 
in the factorization formula (^). For each final state /, we express the de- 
cay amplitudes in terms of parameters a\(f) defined in analogy with similar 
parameters used in the literature on naive factorization. 

4-1 Effective Hamiltonian and decay topologies 
The effective Hamiltonian for B —y Dtt is 

H cS = Q= V* d V ch (CoOo + c s o 8 ) . (15) 

We choose to write the two independent four-quark operators in the singlet- 
octet basis 

O =c7^(1-75)6J 7ai (1-7 5 )u, 

8 = C7"(l - j 5 )T A b d 7M (l - ~f 5 )T A u, (16) 
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rather than in the more conventional basis of 0\ and C>2- The Wilson co- 
efficients Co and C$ describe the exchange of hard gluons with virtualitics 
between the high-energy matching scale Myy and a renormalization scale fj, of 
order mj. (These coefficients are related to the ones of the standard basis by 
Co = C\ + C2/3 and C$ = 2C2.) They are known at next- to- leading order in 
renormalization-group improved perturbation theory and are given by El 

Co = ^t 1 C + + ^^C^, C S = C + -C, (17) 



where 



and 



2N r 2N r 



C „, , = ( 1 + ?M B± \ C±(a») , B ± = ±^p- B , (18) 



~a s (M w )~ 


d± ■ 







a s (M w ) - a s (fi) 
1 + 4^ S± 



(19) 



For N c = 3 and / = 5, we have d + = ^ and d_ = — 1|, as well as S + = 
llll and S- — — . The scheme dependence of the Wilson coefficients at 



next-to-leading order is parameterized by the coefficient B in (J18|). We note 
that -Bndr = 11 in the naive dimensional regularization (NDR) scheme with 
anticommuting 75, and £>hv = 7 in the 't Hooft-Veltman (HV) scheme. We 
will demonstrate below that the scale and scheme dependence of the Wilson 
coefficients is canceled by a corresponding scale and scheme dependence of the 
hadronic matrix elements of the operators Oq and Og. 

Before continuing with a discussion of these matrix elements, it is useful 
to consider the flavour structure for the various contributions to B — > Dn 
decays. The possible quark-level topologies are depicted in Fig. |(| In the 
terminology generally adopted for two-body non-leptonic decays, the decays 
Bd — > D + ir~, Bd — * D°tt° and B~ — > D°ir~ are referred to as class-I, class- 
II and class-Ill, respectively Ej. In Bd — * D + tt~ and B~ — * D°tt~ decays 
the pion can be directly created from the weak current. We call this a class- 

I contribution, following the above terminology. In addition, in the case of 
Bd — > D + tt~ there is a contribution from weak annihilation, and a class- 

II amplitude contributes to B~ — * D°tt~ . The important point is that the 
spectator quark goes into the light meson in the case of the class-II amplitude. 
This amplitude is suppressed in the heavy-quark limit, as is the annihilation 
amplitude. It follows that the amplitude for Bd — > D°ir°, receiving only class- 
II and annihilation contributions, is subleading compared with Bd — > D + tt~ 
and B~ — > D°tt~ , which are dominated by the class-I topology. 
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Figure 6: Basic quark- level topologies for B — * D-k decays (q = u, d): (a) class-I, (b) class- 
II, (c) weak annihilation. Bd — > D + ir~ receives contributions from (a) and (c), Bd — > D°tt° 
from (b) and (c), and B~ — » D°n~ from (a) and (b). Only (a) contributes in the heavy-quark 
limit. 

We shall use the one-loop analysis for B^ — > D + ir~ as a concrete example 
to illustrate explicitly the various steps involved in establishing the factoriza- 
tion formula. Most of the arguments given below are standard from the theory 
of hard exclusive processes involving light hadronsB. However, it is instructive 
to repeat these arguments in the context of B decays. 

4-2 Soft and collinear cancellations at one-loop order 

In order to demonstrate the property of factorization for the decay Bd — ► 
D + tt~, we now analyze the "non-factorizable" one-gluon exchange contribu- 
tions shown in Fig. |5| in some detail. We consider the leading, valence Fock 
state of the emitted pion. This is justified since higher Fock components only 
give power-suppressed contributions to the decay amplitude in the heavy-quark 
limit (as demonstrated later). For the purpose of our discussion, the valence 
Fock state of the pion can be written as 



where a] (b\) denotes the creation operator for a quark (antiquark) in a state 
with spin s =f or s =1, and we have suppressed colour indices. The wave 
function i£(u,l±) is defined as the amplitude for the pion to be composed of 
two on-shell quarks, characterized by longitudinal momentum fraction u and 
transverse momentum lj_. The on-shell momenta of the quark and antiquark 
are chosen as in (g). For the purpose of power counting, l± ~ Aqcd C £~ m^. 
Note that the invariant mass of the valence state is (l q + lq) 2 = l]_/(uu), which 
is of order Aq CD and hence negligible in the heavy-quark limit unless u is 
in the vicinity of the endpoints u = or 1. In this case, the invariant mass 
of the quark-antiquark pair becomes large, and the valence Fock state is no 
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(a\(l q )b\(k)-a\(l q )b\(k)) |0)*(«,il), (20) 
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longer a valid representation of the pion. However, in the heavy-quark limit the 
dominant contributions to the decay amplitude come from configurations where 
both partons are hard (u and u both of order 1), and so the two-particle Fock 
state yields a consistent description. We will provide an explicit consistency 
check of this important feature later on. 

As a next step, we write down the amplitude 

(n(q)\u(0) a d(y) p \0) = f du ^ -J== **(«, f ± ) ( 75 4) a p , (21) 



which appears as an ingredient of the B — > Dir matrix element. It is now 
straightforward to obtain the one-gluon exchange contribution to the B —* Dn 
matrix element of the operator O s . For the sum of the four diagrams in Fig. ^ 
we find 

{D+7r-\0 8 \B d ) 1 . glnon = (22) 



i9 H J (D + \cMk)b\B d ) ¥ ^ _^tr[ 75 4A*{l q ,kM 



where 



7 A (^ C - ft + m c )r r(^ fc + # + m b )Y 



2p c ■ k - fc 2 2p b - k + k 2 



Aa (! 9 , i q , k) = r(fe+ ^ - 2^+ . (23) 

2Z^ • /c -)— fc 2 2/g ■ — )- fc 

Here T = 7^(1 — 75), and p&, p c are the momenta of the b- and c-quark, 
respectively. There is no correction to the matrix element of Oo at order a s , 
because in this case the (du) pair is necessarily in a colour-octet configuration 
and cannot form a pion. 



In (22) the pion wave function ^(u, l±) appears separated from the B — > D 
transition. This is merely a reflection of the fact that we have represented the 
pion state in the form shown in (pp|). It does not, by itself, imply factorization, 
since the right-hand side of (^2|) still involves non-trivial integrations over l± 
and the gluon momentum k, and long- and short-distance contributions are 
not yet disentangled. In order to prove factorization, we need to show that 
the integral over k receives only subdominant contributions from the region 
of small k 2 . This is equivalent to showing that the integral over k does not 
contain infrared divergences at leading power in Aqcd/w;,. 

To demonstrate infrared finitcness of the one-loop integral 

J= jd A k^A 1 (k)®A 2 {l q ,l q ,k) (24) 
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at leading power, the heavy-quark limit and the corresponding large light-cone 
momentum of the pion are again essential. First note that when k is of order 
nib, J ~ 1 by dimensional analysis. Potential infrared divergences could arise 
when k is soft or collinear to the pion momentum q. We need to show that the 
contributions from these regions are power suppressed. (Note that we do not 
need to show that J is infrared finite. It is enough that logarithmic divergences 
have coefficients that are power suppressed.) 

We treat the soft region first. Here all components of k become small 
simultaneously which we describe by scaling k ~ A. Counting powers of A 
(d 4 k ~ A 4 , 1/fc 2 ~ A~ 2 , 1/p-k ~ A -1 ) reveals that each of the four diagrams in 
Fig. corresponding to the four terms in the product in (p4|), is logarithmically 
divergent. However, because k is small the integrand can be simplified. For 
instance, the second term in A2 can be approximated as 



" - If! 4- k\Y 

r, (25) 



2l q ■ k + k 2 2uq-k + 2l±-k+ ^B n - -k + k 2 q ' k 

where we used that to the extreme left or right of an expression gives zero due 
to the on-shcll condition for the external quark lines. We get exactly the same 
expression but with an opposite sign from the other term in A2, and hence the 
soft divergence cancels out. More precisely, we find that the integral is infrared 
finite in the soft region when l± is neglected. When l± is not neglected, there 
is a divergence from soft k which is proportional to l\jm^ ~ Aq CD /to 2 . In 
either case, the soft contribution to J is of order Aqcd/™& or smaller and 
hence suppressed relative to the hard contribution. This corresponds to the 
standard soft cancellation mechanism, which is a technical manifestation of 
colour transparency. 

Each of the four terms in ( p4| ) is also divergent when k becomes collinear 
with the light-cone momentum q. This implies the scaling k + ~ A , k± ~ A, 
and k- - A 2 . Then d 4 fc - dk+ dk~ d 2 'k ± - A 4 , and q ■ k = q+k~ - A 2 , 
k 2 = 2k + k~ + k\ ~ A 2 . The divergence is again logarithmic, and it is thus 
sufficient to consider the leading behaviour in the collinear limit. Writing 
k = aq + . . . we can now simplify the second term of A 2 as 

7a(&+ flr ^ 2(u + a)T _ 
2lq ' k ~\~ k 2 2lq ' k ~\~ k 2 

No simplification occurs in the denominator (in particular, l± cannot be ne- 
glected), but the important point is that the leading contribution is propor- 
tional to q\ . Therefore, substituting k = aq into A\ and using q 2 = 0, we 
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obtain 

. j{j>c + m c )T T(^ b + m b )d 

q\Ai ~ — = , (27) 

2ap c ■ q 2ap b ■ q 

employing the equations of motion for the heavy quarks. Hence the collinear 
divergence cancels by virtue of the standard Ward identity. 

This completes the proof of the absence of infrared divergences at lead- 
ing power in the hard-scattering kernel for Bd — > D + n~ to one-loop order. 
Similar cancellations are observed at higher orders. A complete proof of fac- 
torization at two-loop order can be found ina. Having established that the 
"non-factorizable" diagrams of Fig. ^| are dominated by hard gluon exchange 
(i.e. that the leading contribution to J arises from k of order uib), we may now 
use the fact that <C E to expand A2 in powers of \l±\/E. To leading order 
the expansion simply reduces to neglecting l±_ altogether, which implies l q = uq 
and lq — uq in (^). As a consequence, we may perform the l± integration in 
( p2| ) over the pion distribution amplitude. Defining 

dH ^^ L) -^{u), (28) 



16tt 3 y/2W c 4N C 
the matrix element of Og m (^2|) becomes 

(£»+ 7 r-|0 8 | J B d ) 1 . gluon = (29) 

~- Cy 1 (D+lcAtWblBa) ± U [ du^{u)tT[ l5 4A 2 {uq,uq,k)}- 



SN C J (2tt) 4 x 1 ±v ' 1 0/ k 2 







On the other hand, putting y on the light-cone in (21 ) and comparing with (||), 
we see that the /^-integrated wave function $„•(«) in (^8|) is precisely the light- 
cone distribution amplitude of the pion. This demonstrates the relevance of the 
light-cone wave function to the factorization formula. Note that the collinear 
approximation for the quark and antiquark momenta emerges automatically 
in the heavy-quark limit. 

After the k integral is performed, the expression (|2^) can be cast into the 
form ^ 

(£ + ^|O 8 |^}i- g iuon~^ D (0) / duT s (u,z) $ w (t»), (30) 

Jo 

where z = m c /mb, Tg(u, z) is the hard-scattering kernel, and F B ^ D (0) the 
form factor that parameterizes the (D + |c[. . ]b\Bd) matrix element. Because 
of the absence of soft and collinear infrared divergences in the gluon exchange 
between the (c6) and (du) currents, the hard-scattering kernel T% is calculable 
in QCD perturbation theory. 
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4-3 Matrix elements at next-to-leading order 



We now compute these hard-scattering kernels explicitly to order a s 
effective Hamiltonian (FL5|) can be written as 



The 



G p 
V2 



V* d V c b< 



'Nr. 



1 



2N n 



C+{jm) 



Nr. 



1 



2N r . 



4tt 2N c 



BC s (p) 



O 



(31) 



where the scheme-dependent term in the coefficient of the operator Oq has been 
written explicitly. Because the light-quark pair has to be in a colour singlet 
to produce the pion in the leading Fock state, only Oo gives a contribution 
to zeroth order in a s . Similarly, to first order in a s only Os can contribute. 
The result of evaluating the diagrams in Fig. |^ with an insertion of 0$ can be 
presented in a form that holds simultaneously for a heavy meson H = D, D* 
and a light meson L — 7r, p, using only that the (ud) pair is a colour singlet 
and that the external quarks can be taken on-shell. We obtain (z = m c lrn^) 



(H(p')L(q)\0 8 \B d (p)) 



if L / du§ L {u) 



(32) 



where 



- (61n^ + S) ((J v )-(J A )) + F(u,z) (J v )-F(u,-z)(J A ) 



(Jv) = (H(p')\c(ib\B d (p)), (J A ) = (H{p')\c4 1?> b\B d {p)) . (33) 



It is worth noting that even after computing the one-loop correction the (ud) 
pair retains its V — A structure. This, together with (^), implies that the form 
of (^) is identical for pions and longitudinally polarized p mesons. (The pro- 
duction of transversely polarized p mesons is power suppressed in Aqcd /mb-) 
The function F(u, z) appearing in (|32]) is given by 



where 
f(u,z) 



F(u, z) = (z + 2 In |) In z 2 - 7 + f(u, z) + f(u, l/z) , 
u(l - z 2 )[3(l - u(l - z 2 )) + z] 



(34) 



[l-u(l- 



-2M2 



\n[u(l - z 2 )] 
l-u(l-z 2 ) 



- ln 2 [u(l - z 2 )] - Li 2 [l - u(l - z 2 )] - {i 



(35) 
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and Li2(x) is the dilogarithm. The contribution of f(u,z) in (34) comes from 
the first two diagrams in Fig. |B| with the gluon coupling to the b quark, whereas 
f(u,l/z) arises from the last two diagrams with the gluon coupling to the 
charm quark. Note that the terms in the large square brackets in the definition 
of the function f(u, z) vanish for a symmetric light-cone distribution amplitude. 
These terms can be dropped if the light final-state meson is a pion or a p 
meson, but they are relevant, e.g., for the discussion of Cabibbo-suppressed 
decays such as B d -> D^+K~ and B d D^+K*~. 

The discontinuity of the amplitude, which is responsible for the occurrence 
of the strong rescattering phase, arises from f(u, 1/z) and can be obtained by 
recalling that z 2 is z 2 — it with e > infinitesimal. We find 



• Im F(u, z) 



(1 - tt)(l - z 2 )[3(l - u(l - z 2 )) + z] 



[l-u(l-z 2 )} 2 



ln[l-it(l-z 2 )] + 21nw + 



l-u(l-z 2 ) 



{u — > it} 



(36) 



As mentioned above, (|3^) is applicable to all decays of the type B d — ^ 
D^+L~, where L is a light hadron such as a pion or a (longitudinally polar- 
ized) p meson. Only the operator Jy contributes to Bd — » D + L~, and only 
J A contributes to Bd — > D* + L~ . Our result can therefore be written as 



(D + L-\O , 8 \Bd) = (D+\c^(l- l5 )b\B d ) -if L q^ duT 0i8 (u,z)^ L (u), (37) 

Jo 

where L = n, p, and the hard-scattering kernels are 
T (u,z) = l + 0{a 2 s ), 



T 8 (u, z) 



47 2N~ r 



-6 In 



B + F(u,z) 



+ 0(a 2 ) . 



(38) 



When the D meson is replaced by a D* meson, the result is identical except 
that F(u,z) must be replaced with F(u, — z). Since no order-a s corrections 
exist for Oq, the matrix element retains its leading-order factorized form 



(D + L-\O \B d ) = if Lqil (D+\c^(l - j 5 )b\B d ) 



(39) 



to this accuracy. From ( pq ) it follows that Tg(u, z) tends to a constant as u ap- 
proaches the endpoints (u — > 0, 1). (This is strictly true for the part of T 8 (u, z) 
that is symmetric in u «-> u; the asymmetric part diverges logarithmically as 
u — > 0, which however does not affect the power behaviour and the convergence 
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properties in the endpoint region.) Therefore the contribution to ( p7| ) from the 
endpoint region is suppressed, both by phase space and by the endpoint sup- 
pression intrinsic to $i(it). Consequently, the emitted light meson is indeed 
dominated by energetic constituents, as required for the self-consistency of the 
factorization formula. 

The final result for the class-I, non-leptonic Bd — > D" + L~ decay am- 
plitudes, in the heavy-quark limit and at next-to- leading order in a s , can be 
compactly expressed in terms of the matrix elements of a "transition operator" 



T 



V2 



v: d v cb 



a\(DL) Q v - ai (D*L) Q A 



(40) 



where 

Qv = cY L b <g> d7 M (l - 75)" ! 



Qa =07^756 <8> rf7 M (l-7 5 )u, (41) 



and hadronic matrix elements of Qy.A are understood to be evaluated in fac- 
torized form, i.e. 



{DL\j x ®n\B) = {D\j x \B) (L|j 2 |0> 



(42) 



Eq. (|ic| ) defines the quantities ai(D^L), which include the leading "non- 
factorizable" corrections, in a renormalization-scale and -scheme independent 
way. To leading power in Aqcd /nib these quantities should not be interpreted 
as phenomenological parameters (as is usually done), because they are domi- 
nated by hard gluon exchange and thus calculable in QCD. At next-to-leading 
order we get 



a x {DL) 



N c + \ 



2N L 



N,. 



2N n 



-C-(jj) 



a x {D*L) 



Att2N< 



c 8 (m) 



61nAj+/ duF(u,z)<f> L (u) 
m t Jo 



2N L 

a s C F 
4tt 2N { 



C+(p) 



Nr. 



2N n 



-C-(jj) 



c 8 (/0 



-61n-^+ / duF(u,-z)<f> L (u) 
mt Jo 



(43) 



We observe that the scheme-dependent terms parameterized by B have can- 
celed between the coefficient of Oq in ( |3l| ) and the matrix element of 0$ in 
([57]). Likewise, the fj, dependence of the terms in brackets in (|I^) cancels 
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against the scale dependence of the coefficients C±(/x), ensuring a consistent 
result at next-to-leading order. The coefficients ai(DL) and a\(D*L) are seen 
to be non-universal, i.e. they depend explicitly on the nature of the final-state 
mesons. This dependence enters via the light-cone distribution amplitude of 
the light emission meson and via the analytic form of the hard-scattering kernel 
(F(u,z) vs. F(u,—z)). However, the non-universality enters only at next-to- 
leading order. 

Using the fact that violations of heavy-quark spin symmetry require hard 
gluon exchange, Politzer and Wise have computed the "non-factorizable" ver- 
tex corrections to the decay-rate ratio of the Dir and D*ir final states many 
years ago 113. In the context of our formalism, this calculation requires the 
symmetric part (with respect to u <-> u) of the difference F(u,z) — F(u,—z). 
Explicitly, 



T(B d -> D+7T-) 
T{B d -> D*+tt-) 



(D+\c<i(l- l5 )b\B d ) 



(D*+\c4(l- lh )b\B d ) 



ai(Dir) 



ai(Z)*7r) 



(44) 



where for simplicity we neglect the light meson masses as well as the mass 
difference between D and D* in the phase-space for the two decays. At next- 
to-leading order 



ai{Dir) 



ai(£>*7r) 



1 + ?^fSJ Re / du[F{u,z)-F{u,-z)]^{u). (45) 
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Our result for the symmetric part of F(u, z) — F(u, — z) coincides with that 
found in 113. 



5 Power-suppressed contributions 

Up to this point we have presented arguments in favour of factorization of 
non-leptonic i?-decay amplitudes in the heavy-quark limit, and have explored 
in detail how the factorization formula works at one-loop order for the decays 
B d -► DM+L-. It is now time to show that other contributions not considered 
so far are indeed power suppressed. This is necessary to fully establish the 
factorization formula. Besides, it will also provide some numerical estimates 
of the corrections to the heavy-quark limit. 

We start by discussing interactions involving the spectator quark and weak 
annihilation contributions, before turning to the more delicate question of the 
importance of non- valence Fock states. 
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Figure 7: "Non-factorizable" spectator interactions. 



5.1 Interactions with the spectator quark 

Clearly, the diagrams shown in Fig. fj] cannot be associated with the form- 
factor term in the factorization formula (||). We will now show that for B 
decays into a heavy-light final state their contribution is power suppressed in 
the heavy-quark limit. (This suppression does not occur for decays into two 
light mesons, where hard spectator interactions contribute at leading power. 
In this case, they contribute to the kernels T/ in the factorization formula 
(second term in Fig. [j]).) 

In general, "non-factorizable" diagrams involving an interaction with the 
spectator quark would impede factorization if there existed a soft contribution 
at leading power. While such terms are present in each of the two diagrams 
separately, they cancel in the sum over the two gluon attachments to the (ud) 
pair by virtue of the same colour-transparency argument that was applied to 
the "non-factorizable" vertex corrections. 

Focusing again on decays into a heavy and a light meson, such as Bd ~^ 
D + n~, we still need to show that the contribution remaining after the soft 
cancellation is power suppressed relative to the leading-order contribution ([l3]) . 
A straightforward calculation leads to the following (simplified) result for the 
sum of the two diagrams: 



A(B d 



D+ir- 



spec ~ Gf fnfDfB «s 



io ? Jo V Jo u 

~ G F a s m b Aq CD . (46) 

This is indeed power suppressed relative to (|l3|) . Note that the gluon virtuality 
is of order ^wrn^ ~ Aq CD and so, strictly speaking, the calculation in terms 
of light-cone distribution amplitudes cannot be justified. Nevertheless, we use 
( [46| ) to deduce the scaling behaviour of the soft contribution, as we did for the 
heavy-light form factor in Sect. 3.5. 



28 



^ >4 

(a) (b) (c) (d) 

Figure 8: Annihilation diagrams. 

5.2 Annihilation topologies 

Our next concern are the annihilation diagrams shown in Fig. || which also 
contribute to the decay Bd — ► D + ir~ . The hard part of these diagrams could, 
in principle, be absorbed into hard-scattering kernels of the type T/ 7 . The soft 
part, if unsuppressed, would violate factorization. However, we will see that 
the hard part as well as the soft part are suppressed by at least one power of 
Aqcd/to6- 

The argument goes as follows. We write the annihilation amplitude as 

A{Bd — * -D + 7r~)aim — Gf fwfofB a s 

x [ d^drjdu $s (£) $d (v) (u) T" uin (£, 77, u) , (47) 
Jo 

where the dimensionless function T ann (£, 77, u) is a product of propagators and 
vertices. The product of decay constants scales as Aq CD /mb. Since d£ $b (£ ) 
scales as 1 and so does dij $15(77), while du$> n (u) is never larger than 1, the 
amplitude can only compete with the leading-order result ( |l3| ) if T ann (£, 77, u) 
can be made of order (mj/ Aqcd) 3 or larger. Since T ann (£, 77, it) contains only 
two propagators, this can be achieved only if both quarks the gluon splits into 
are soft, in which case T ann (£, 77, u) ~ (™&/Aqcd) 4 - But then dw$ T (u) ~ 
(Aqcd/to6) 2 , so that this contribution is power suppressed. 

5.3 Non-leading Fock states 

Our discussion so far concentrated on contributions related to the quark 
antiquark components of the meson wave functions. We now present quali- 
tative arguments that justify this restriction tOj-the valence-quark Fock com- 
ponents. Some of these arguments are standard Era. We will argue that higher 
Fock states yield only subleading contributions in the heavy-quark limit. 
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Figure 9: Diagram that contributes to the hard-scattering kernel involving a quark- 
antiquark-gluon distribution amplitude of the B meson and the emitted light meson. 



Additional hard partons 

An example of a diagram that would contribute to a hard-scattering function 
involving quark-antiquark-gluon components of the emitted meson and the 
B meson is shown in Fig. ^. For light mesons, higher Fock components are 
related to higher-order terms in the collinear expansion, including the effects 
of intrinsic transverse momentum and off-shcllness of the partons by gauge 
invariance. The assumption is that the additional partons are collinear and 
carry a finite fraction of the meson momentum in the heavy-quark limit. Under 
this assumption, it is easy to see that adding additional partons to the Fock 
state increases the number of off-shell propagators in a given diagram (compare 
Fig. [)]to Fig. H). This implies power suppression in the heavy-quark expansion. 
Additional partons in the _B-meson wave function are always soft, as is the 
spectator quark. Nevertheless, when these partons are connected to the hard- 
scattering amplitudes the virtuality of the additional propagators is still of 
order m^AQCD, which is sufficient to guarantee power suppression. 




Figure 10: The contribution of the qqg Fock state to the —> D + tt decay amplitude. 



Let us study in more detail how the power suppression arises for the sim- 
plest non-trivial example, where the pion is composed of a quark, an antiquark, 
and an additional gluon. The contribution of this 3-particle Fock state to the 
B — > Dtt decay amplitude is shown in Fig. 10. It is convenient to use the 



Fock-Schwinger gauge, which allows us to express the gluon field A\ in terms 
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of the field-strength tensor G p \ via 



A x (x)= [ dvvx p G pX {vx). (48) 
Jo 

Up to twist-4 level, there are three quark-antiquark-gluon matrix elements 
that could potentially contribute to the diagrams shown in Fig. [l0| Due to the 
V — A structure of the weak-interaction vertex, the oplv-jjelevant three-particle 
light-cone wave function has twist-4 and is given byLLaH 

(n(q)\d(0) 7 ^ 5 g s G al 3(vx)u(0)\0) 

+ U — (Qaxp - q x a ) f Vu (M^) + H e* "**-* . (49) 
q-x J 

Here J T>u = f Q dui du 2 <f«3(5(l — Mi — — u 3 ), with u\, u 2 and m 3 the fractions 
of the pion momentum carried by the quark, antiquark and gluon, respectively. 
Evaluating the diagrams in Fig. |l^, and neglecting the charm-quark mass for 
simplicity, we find 

(D + TT-\O s \B d ) qqg =iU{D+\ci(l-^)b\B d ) [ Vu 2 ^ U ^ . (50) 

Since <fi\\ ~ Aq CD , the suppression by two powers of AQcn/mb compared to 
the leading-order matrix element is obvious. Note that due to G-parity <f>u is 
antisymmetric in u\ <-*■ u 2 for a pion, so that (pfj) vanishes in this case. 



Additional soft partons 

A more precarious situation may arise when the additional Fock components 
carry only a small fraction of the mesommomentum, contrary to the assump- 
tion made above. It is usually argued Em that these configurations are sup- 
pressed, because they occupy only a small fraction of the available phase space 
(since J dui ~ Aqcd / mb when the parton that carries momentum fraction u, is 
soft). This argument does not apply when the process involves heavy mesons. 
Consider, for example, the diagram shown in Fig. ^ (a) for the decay B — > Dir. 
Its contribution involves the overlap of the B-meson wave function involving 
additional soft gluons with the wave function of the D meson, also containing 
soft gluons. There is no reason to suppose that this overlap is suppressed rela- 
tive to the soft overlap of the valence-quark wave functions. It represents (part 
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Figure 11: (a) Soft overlap contribution which is part of the B — » D form factor, (b) Soft 
overlap with the pion which would violate factorization, if it were unsuppressed. 

of) the overlap of the "soft cloud" around the b quark with (part of) the "soft 
cloud" around the c quark after the weak decay. The partonic decomposition 
of this cloud is unrestricted up to global quantum numbers. (In the case where 
the B meson decays into two light mesons, there is a form-factor suppression 
~ (AQCD/wib) 3 / 2 for the overlap of the valence-quark wave functions, but once 
this price is paid there is again no reason for further suppression of additional 
soft gluons in the overlap of the i?-meson wave function and the wave function 
of the recoiling meson.) 

The previous paragraph essentially repeated our earlier argument against 
the hard-scattering approach, and in favour of using the B —> D form factor as 
an input to the factorization formula. However, given the presence of additional 
soft partons in the B — > D transition, we must now argue that it is unlikely 
that the emitted pion drags with it one of these soft partons, for instance a soft 
gluon that goes into the pion wave function, as shown in Fig. [ll] (b). Notice 
that if the (qq) pair is produced in a colour-octet state, at least one gluon 
(or a further (qq) pair) must be pulled into the emitted meson if the decay 
is to result in a two-body final state. What suppresses the process shown in 
Fig. [ll] (b) relative to the one in Fig. [ll] (a) even if the emitted (qq) pair is in 
a colour-octet state? 

It is once more colour transparency that saves us. The dominant config- 
uration has both quarks carry a large fraction of the pion momentum, and 
only the gluon might be soft. In this situation we can apply a non-local "op- 
erator product expansion" to determine the coupling of the soft gluon to the 
small (qq) paira. The gluon endpoint behaviour of the qqg wave function is 
then determined by the sum of the two diagrams shown on the right-hand 
side in Fig. ^[ The leading term (for small gluon momentum) cancels in the 
sum of the two diagrams, because the meson (represented by the black bar) 
is a colour singlet. This cancellation, which is exactly the same cancellation 
needed to demonstrate that "non-factorizable" vertex corrections are domi- 
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Figure 12: Quark— antiquark-gluon distribution amplitude in the gluon endpoint region. 



nated by hard gluons, provides one factor of Aqqd /rrib needed to show that 
Fig. [ll] (b) is power suppressed relative to Fig. [n] (a). 

In summary, we have (qualitatively) covered all possibilities for non- valence 
contributions to the decay amplitude and find that they are all power sup- 
pressed in the heavy-quark limit. 

6 Limitations of the factorization approach 

The factorization formula (^) holds in the heavy-quark limit m& — > oo. Cor- 
rections to the asymptotic limit are power-suppressed in the ratio Aqcd/w;, 
and, generally speaking, do not assume a factorized form. Since is fixed to 
about 5 GeV in the real world, one may worry about the magnitude of power 
corrections to hadronic -B-decay amplitudes. Naive dimensional analysis would 
suggest that these corrections should be of order 10% or so. We now discuss 
several reasons why some power corrections could turn out to be numerically 
larger than suggested by the parametric suppression factor Aqcb/itt-i,. Most 
of these "dangerous" corrections occur in more complicated, rare hadronic B 
decays into two light mesons, but are absent in decays such as B — > Dir. 

6. 1 Several small parameters 

Large non-factorizable power corrections may arise if the leading-power, fac- 
torizable term is somehow suppressed. There are several possibilities for such 
a suppression, given a variety of small parameters that may enter into the 
non-leptonic decay amplitudes. 

i) The hard, "non-factorizable" effects computed using the factorization for- 
mula occur at order a s . Some other interesting effects such as final-state 
interactions appear first at this order. For instance, strong-interaction 
phases due to hard interactions are of order a s , while soft rescattering 
phases are of order Aqcd/w&- Since for realistic B mesons a s is not par- 
ticularly large compared to Aqcd/to&, we should not expect that these 
phases can be calculated with great precision. In practice, however, it is 
probably more important to know that the strong-interaction phases are 
parametrically suppressed in the heavy-quark limit and thus should be 
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small. (This does not apply if the real part of the decay amplitudes is 
suppressed for some reason; see below.) 

ii) If the leading, lowest-order (in a s ) contribution to the decay amplitude 
is colour suppressed, as occurs for the class-II decay Ba ~ * 7r°7r°, then 
perturbative and power corrections can be sizeable. In such a case even 
the hard strong-interaction phase of the amplitude can be largeEra. But 
at the same time soft contributions could be potentially important, so 
that in some cases only an order-of-magnitude estimate of the amplitude 
may be possible. 

iii) The effective Hamiltonian ([!]) contains many Wilson coefficients C, that 
are small relative to C\ ~ 1. There are decays for which the entire 
leading-power contribution is suppressed by small Wilson coefficients, 
but some power-suppressed effects are not. An example of this type is 
B~ — ► K~K°. This decay proceeds through a penguin operator b — » dss 
at leading power. But the annihilation contribution, which is power 
suppressed, can occur through the current-current operator with large 
Wilson coefficient C\. Our approach does not apply to such (presum- 
ably) annihilation-dominated decays, unless a systematic treatment of 
annihilation amplitudes can be found. 

iv) Some amplitudes may be suppressed by a combination of small CKM ma- 
trix elements. For example, B — > nK decays receive large penguin con- 
tributions despite their small Wilson coefficients, because the so-called 
tree amplitude is CKM suppressed. This is not a problem for factoriza- 
tion, since it applies to the penguin and the tree amplitudes. We are 
not aware of any case (for ordinary B mesons) in which a purely power- 
suppressed term is CKM enhanced and which would therefore dominate 
the decay amplitude. (But this situation could occur for B~ — > D°K~, 
where the QCD dynamics is similar if we consider the charm quark as a 
light quark.) 



6.2 Power corrections enhanced by small quark masses 

There is another enhancement of power-suppressed effects for some decays into 
two light mesons, connected with the curious numerical fact that 

2ml Man) , . 
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is much larger than its naive scaling estimate Aqcd- (Here (qq) = (0\uu\0) = 
(0|dd|0) is the quark condensate.) Consider the contribution of the penguin 



34 



operator Oq — (dibj)v-A(ujUi)v+A to the Bd — > -k + tt decay amplitude. The 
leading-order graph of Fig. || results in the expression 



(ir + TT~\(dibj)v—A{ujUi)v+A\Bd) = im 2 B F^(0)U x — , (52) 

which is formally a Aqcd / mb power correction compared to the corresponding 
matrix element of a product of two left-handed currents, but numerically large 
due to ([si]) . We would not have to worry about such terms if they could all 
be identified and the factorization formula (^) applied to them, since in this 
case higher-order perturbative corrections would not contain non-factorizing 
infrared logarithms. However, this is not the case. 

After including radiative corrections, the matrix element on the left-hand 
side of ([52]) is expressed as a non-trivial convolution with pion light-cone dis- 
tribution amplitudes. The terms involving \i % can be related to two-particle 
twist-3 (rather than leading twist-2) distribution amplitudes, conventionally 
called and $ CT (it). We find that the radiative corrections to the matrix 

element in (^2|) do indeed factorize. However, at the same order there ap- 
pear twist-3 corrections to the hard spectator interaction shown in Fig. 0, and 
these contributions contain an endpoint divergence (related to the fact that the 
distribution amplitudes & p (u) and &' a (u) do not vanish at the endpoints). In 
other words, the twist-3 "corrections" to the hard spectator term in the second 
factorization formula in (^) relative to the "leading" twist-2 contributions are 
of the form a s x logarithmic divergence, which we interpret as being of order 
1. The non-factorizing character of the "chirally-enhanced" power corrections 
can introduce a substantial uncertainty ia. some decay modes 113. As in the 
related situation for the pion form factored, one may argue that the endpoint 
divergence is suppressed by a Sudakov form factor. However, it is likely that 
when 7ti 6 is not large enough to suppress these chirally-enhanced terms, then 
it is also not large enough to make Sudakov suppression effective. 

We stress that the chirally-enhanced terms do not appear in decays into 
a heavy and a light meson such as B — > Dtt, because these decays have no 
penguin contribution and no contribution from the hard spectator interaction. 
Hence, the twist-3 light-cone distribution amplitudes responsible for chirally- 
enhanced power corrections do not enter in the evaluation of the decay ampli- 
tudes. 

6.3 Non-leptonic decays when M 2 is not light 



The analysis of non-leptonic decay amplitudes in Sect. 3.3 referred to decays 
where the emission particle M2 is a light meson. We now briefly discuss the 
case where M 2 is heavy. 
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Suppose that M2 is a D meson, whereas the meson that picks up the 
spectator quark can be heavy or light. Examples of this type are the decays 
Bd — > tt°D° and Bd — * D + D~ . It is intuitively clear that factorization must 
be problematic in these cases, because the heavy D meson has a large overlap 
with the Bit or BD systems, which are dominated by soft processes. In more 
detail, we consider the coupling of a gluon to the two quarks that form the 
emitted D meson, i.e. the pairs of diagrams in Fig. I (a+b), (c+d) and Fig. @. 
Denoting the gluon momentum by k, the quark momenta by l q and l q , and the 
D-meson momentum by q, we find that the gluon couples to the "current" 

^ = lx{y q + fj + m c )T _ T(f/ q + #) 7A 
2Zq * k ~\~ kP 2Zg ■ k ~~\~ kP" 

where T is part of the weak decay vertex. When k is soft (all components of 
order Aqcd), each of the two terms scales as 1/Aqqd- Taking into account 



the complete amplitude as done explicitly in Sect. 4.2, we can see that the 
decoupling of soft gluons requires that the two terms in ( |53"| ) cancel, leaving a 
remainder of order l/m&. This cancellation does indeed occur when M2 is a 
light meson, since in this case l q and l q are dominated by their longitudinal 
components. When M2 is heavy, the momenta l q and lq are asymmetric, with 
all components of the light antiquark momentum l q of order Aqcd in the B- 
or D- meson rest frames, while the zero-component of l q is of order mj,. Hence 
the current can be approximated by 

k 21 q ■ k + k z Aqcd 

and the soft cancellation does not occur. (The on-shell condition for the charm 
quark has been used to arrive at this equation.) 

It follows that the emitted D meson does not factorize from the rest of the 
process, and that a factorization formula analogous to (||) does not apply to 
decays such as Bd — * tt°D° and Bd — > D + D~ . An important implication of 
this statement is that one should also not expect naive factorization to work 
in these cases. In other words, we expect that "non-factorizable" corrections 
modify the factorized decay amplitudes by terms of order 1. 

6.4 Difficulties with charm 

There are decay modes, such as B~ — > D°ir~ , in which the spectator quark 
can go to either of the two final-state mesons. The factorization formula ([|) 
applies to the contribution that arises when the spectator quark goes to the D 
meson, but not when the spectator quark goes to the pion. However, even in 
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the latter case we may use naive factorization to estimate the power behaviour 
of the decay amplitude. Adapting ( |i"3"| ) to the decay B~ — ► D°tt~ , we find 
that the non-factorizing (class-II) amplitude is suppressed compared to the 
factorizing (class-I) amplitude by 

A(B- ^ZJ°7r-) class -„ F B ^(m 2 D )f D / Aqqd \ 2 

Here we use that F B ^ 7r (q 2 ) ~ (Aqcd/w;>) 3 / 2 even for g 2 ~ to 2 ,, as long as 
*?ma£c — 9 2 i s a l so °f order to 2 . (It follows from our definition of heavy final-state 
mesons that these conditions are fulfilled.) As a consequence, strictly speaking 
factorization does hold for B~ — > D°ir~ decays in the sense that the class-II 
contribution is power suppressed with respect to the class-I contribution. 

Unfortunately, the scaling behaviour for real B and D mesons is far from 
the estimate (^H|) valid in the heavy-quark limit. Based on the dominance of 
the class-I amplitude we would expect that 

Br(B d -> D+TT-) V ; 

in the heavy-quark limit. This contradicts existing data which yield R = 
1.89 ±0.35, despite the additional colour suppression of the class-II amplitude. 
One reason for the failure of power counting lies in the departure of the decay 
constants and form factors from naive power counting. The following compares 
the power counting to the actual numbers (square brackets): 

/d /AgoD^' 2 r ^ ^|)^^ a/2 N0 . 5] . (57) 



U \ m c J Ff^ D (ml) V m b 

However, it is unclear whether the failure of power counting can be attributed 
to the form factors and decay constants alone. 

Note that for the purposes of power counting we treated the charm quark 
as heavy, taking the heavy-quark limit for fixed m c /mb. This simplified the 
discussion, since we did not have to introduce m c as a separate scale. However, 
in reality charm is somewhat intermediate between a heavy and a light quark, 
since m c is not particularly large compared to Aqcd- In this context it is 
worth noting that the first hard-scattering kernel in (^) cannot have Aqcd/™c 
corrections, since there is a smooth transition to the case of two light mesons. 
The situation is different with the hard spectator interaction term, which we 
argued to be power suppressed for decays into a D meson and a light meson. 



We shall come back to this in Sect. 7.5, where we estimate the magnitude of 



this term for the Dn final state, relaxing the assumption that the D meson is 
heavy. 
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7 Phenomenology of B — ► £)(*)_£, decays 



The matrix elements we have computed in Sect. 4.3 provide the theoretical 
basis for a model-independent calculation of the class-I non-leptonic decay 
amplitudes for decays of the type B — ► D^L, where L is a light meson, to 
leading power in Aqcd/™6 and at next-to-leading order in renormalization- 
group improved perturbation theory. In this section we discuss phenomeno- 
logical applications of this formalism and confront our numerical results with 
experiment. We also provide some numerical estimates of power-suppressed 
corrections to the factorization formula. 

7.1 Non-leptonic decay amplitudes 

The results for the class-I decay amplitudes for B — > D^L are obtained by 
evaluating the (factorized) hadronic matrix elements of the transition opera- 
tor T defined in ( phf ) . They are written in terms of products of CKM matrix 
elements, light-meson decay constants, B — ► transition form factors, and 
the QCD parameters a\{D^L). The decay constants can be determined ex- 
perimentally using data on the weak leptonic decays P~ — ► l~Pi( , y), hadronic 
t~ j=f M~v r decays, and the electromagnetic decays V° — > e + e~. Follow- 
inglia, we use U = 131 MeV, f K = 160 MeV, f p = 210 McV, f K * = 214 MeV, 
and f ai = 229 MeV. (Here a\ is the pseudovector meson with mass m ai ~ 
1230 MeV.) 

The non-leptonic Bd — > D^ + L~ decay amplitudes for L — n, p can be 
expressed as 

A(B d -» D+n-) = v: d V cb ai (Dn) U F (m*) (m| - m%) , 
A(B d ^D*+ir-) = -i^V: d V cbai (D*n)UA (ml)2m D .s*-p, 



A(B d D+p-) = -«% V* d Vcb Oi (Dp) f p F+{m 2 p ) 2m p rf ■ p , 



(58) 



where p (p') is the momentum of the B (charm) meson, s and rj are polarization 
vectors, and the form factors Fo, F + and A$ are defined in the usual waytLS. 
The decay mode B d — > D* + p~ has a richer structure than the decays with at 
least one pseudoscalar in the final state. The most general Lorentz-invariant 
decomposition of the corresponding decay amplitude can be written as 

A(B d -> D*+p-) = & v; d V cb e*»rf» (Sx g pv - S 2 q p p' u + iS 3 e^ a/3 p'Vj , 

(59) 
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where the quantities Si can be expressed in terms of semi-leptonic form factors. 
To leading power in Aqcd/wj,, we obtain 

5 1 = a 1 (D*p)m p f p (m B + m D *)A 1 (m 2 ) , 

2A 2 (m 2 ) 

5 2 = ai (D*p) m p f p 2 -±-£>- . (60) 

m B + tti d* 

The contribution proportional to S3 in ( p9| ) is associated with transversely 
polarized p mesons and thus leads to power-suppressed effects, which we do 
not consider here. 

The various B — > form factors entering the expressions for the de- 
cay amplitudes can be determined by combining experimental data on semi- 
leptonic .cLecays with theoretical relations derived using heavy-quark effective 
theory HO. Since we work to leading order in Aqcd/^6, it is consistent to set 
the light meson masses to zero and evaluate these form factors at q 2 = 0. In 
this case the kinematic relations 

F (0) = F+(0) , (m B + m D .)Ai(0) - (m B - rn D *)A 2 (0) = 2m D ,A Q (0) 

(61) 

allow us to express the two Bd — > D + L rates in terms of F+(0), and the two 
Bd — > D* + L~ rates in terms of Aq(0). Heavy-quark symmetry implies that 
these two form factors are equal to within a few percent til. Below we adopt 
the common value F+(0) = ^4o(0) = 0.6. All our predictions for decay rates 
will be proportional to the square of this number. 

7.2 Meson distribution amplitudes and predictions for a\ 

Let us now discuss in more detail the ingredients required for the numerical 
analysis of the coefficients a\{D^*>L). The Wilson coefficients C, in the effective 
weak Hamiltonian depend on the choice of the scale p as well as on the value 
of the strong coupling a s , for which we take ct s (m,z) — 0.118 and two-loop 
evolution down to a scale p ~ mj,. To study the residual scale dependence 
of the results, which remains because the perturbation series are truncated at 
next-to-leading order, we vary p between mb/2 and 2mf,. The hard-scattering 
kernels depend on the ratio of the heavy-quark masses, for which we take 
z = m c /m b = 0.30 ±0.05. 

Hadronic uncertainties enter the analysis also through the parameteriza- 
tions used for the meson light-cone distribution amplitudes. It is convenient 
and conventional to expand the distribution amplitudes in Gegenbauer poly- 
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Table 1: Numerical values for the integrals J du F(u, z) (u) (upper portion) and 
1 duF(u, —z) (lower portion) obtained including the first two Gegenbauer moments. 



z 


Leading term 


Coefficient of a\ 


Coefficient of a\ 


0.25 
0.30 
0.35 


-8.41-9.51i 
-8.79- 9.09i 
-9.13-8.59i 


5.92- 12.19i 
5.78- 12. 7H 
5.60- 13.21* 


-1.33 + 0.36t 
-1.19 + 0.58i 
-1.00 + 0.73t 


0.25 
0.30 
0.35 


-8.45- 6.56i 
-8.37- 5. 99i 
-8.24- 5.44i 


6.72- 10.73* 
6.83- 11.49i 
6.81 - 12.29i 


-0.38 + 0.93i 
-0.21 + 0.85i 
-0.08 + 0.75i 



nomials as 



= 6u(l - u) 



(62) 



where C[ 3 ^ 2 \x) = 3x, C% 3 ' 2 '(x) — |(5a; 2 — 1), etc. The Gegenbauer mo- 
ments o^(/i) are multiplicatively renormalized. The scale dependence of these 
quantities would, however, enter the results for the coefficients only at order 
a 2 , which is beyond the accuracy of our calculation. We assume that the 
leading-twist distribution amplitudes are close to their asymptotic form and 
thus truncate the expansion at n = 2. However, it would be straightforward 
to account for higher-order terms if desired. For the asymptotic form of the 
distribution amplitude, <&l(u) = 6tt(l — u), the integral in M3) yields 



duF(u,z)$ L (u) = 31nz 2 -7 



6z(l-2z) 
(l-z) 2 (l + z) 3 



Li 2 (z 2 ) 



3(2-3z + 2z 2 + z 3 ) 



{l-z){l + zf 



4- 17z + 20z 2 + 5z 3 
2(l-z)(l + z) 2 



+ {z^l/z} 



ln(l - z 2 ) 



(63) 



and the corresponding result with the function F(u, —z) is obtained by re- 
placing z — > — z. More generally, a numerical integration with a distribution 
amplitude expanded in Gegenbauer polynomials yields the results collected in 
Table [t| We observe that the first two Gegenbauer polynomials in the ex- 
pansion of the light-cone distribution amplitudes give contributions of similar 
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Table 2: The QCD coefficients ai(£H*) L) at next-to-leading order for three different values 
of the renormalization scale fi. The leading-order values are shown for comparison. 





\i = m b /2 


H = m b 


fi = 2m;, 


ai(DL) 


1.074 + 0.037* 


1.055 + 0.020? 


1.038 + 0.011? 




-(0.024- 0.052i)af 


-(0.013 -0.028i) of 1 


-(0.007- 0.015?) 


ai{D*L) 


1.072 + 0.024i 


1.054 + 0.013? 


1.037 + 0.007? 




-(0.028- 0.047i)af 


-(0.015- 0.025?) 


-(0.008- 0.014?) af 


a\° 


1.049 


1.025 


1.011 



magnitude, whereas the second moment gives rise to much smaller effects. This 
tendency persists in higher orders. For our numerical discussion it is a safe 
approximation to truncate the expansion after the first non-trivial moment. 
The dependence of the results on the value of the quark mass ratio z = m c /mi, 
is mild and can be neglected for all practical purposes. We also note that the 
difference of the convolutions with the kernels for a pseudoscalar D and vector 
D* meson are numerically very small. This observation is, however, specific to 
the case of B — > D^L decays and should not be generalized to other decays. 

Next we evaluate the complete results for the parameters a\ at next-to- 
lcading order, and to leading power in Aqcd/w&. We set z — m c jm b = 0.3. 
Varying z between 0.25 and 0.35 would change the results by less than 0.5%. 
The results are shown in Table |[ The contributions proportional to the second 
Gegenbauer moment a\ have coefficients of order 0.2% or less and can safely 
be neglected. The contributions associated with a\ are present only for the 
strange mesons K and K*, but not for n and p. Moreover, the imaginary 
parts of the coefficients contribute to their modulus only at order a 2 s , which is 
beyond the accuracy of our analysis. To summarize, we thus obtain 

\a x {DL)\ = 1.055t° ^ 9 7 - (0.013l^K » 

\ ai (D*L)\ = 1.054ir o i 8 7 - (0.015±roo?K , (64) 

where the quoted errors reflect the perturbative uncertainty due to the scale 
ambiguity (and the negligible dependence on the value of the ratio of quark 
masses and higher Gegenbauer moments), but not the effects of power-sup- 
pressed corrections. These will be estimated later. It is evident that within 
theoretical uncertainties there is no significant difference between the two a\ 
parameters, and there is only a very small sensitivity to the differences between 
strange and non-strange mesons (assuming that ) \ < 1). In our numerical 
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analysis below we thus take |<xi| = 1.05 for all decay modes. 



7.3 Tests of factorization 

The main lesson from the previous discussion is that corrections to naive fac- 
torization in the class-I decays B d — ► D^ + L~ are very small. The reason 
is that these effects are governed by a small Wilson coefficient and, more- 
over, are colour suppressed by a factor 1/N 2 . For these decays, the most 
important implications of the QCD factorization formula are to restore the 
renormalization-group invariance of the theoretical predictions, and to pro- 
vide a theoretical justification for why naive factorization works so well. On 
the other hand, given the theoretical uncertainties arising, e.g., from unknown 
power-suppressed corrections, there is little hope to confront the extremely 
small predictions for non-universal (process-dependent) "non-factorizable" cor- 
rections with experimental data. Rather, what we may do is ask whether data 
supports the prediction of a quasi- universal parameter |eti| ~ 1.05 in these 
decays. If this is indeed the case, it would support the usefulness of the heavy- 
quark limit in analyzing non-leptonic decay amplitudes. If, on the other hand, 
we were to find large non-universal effects, this would point towards the exis- 
tence of sizeable power corrections to our predictions. We will see that within 
present experimental errors the data are in good agreement with our prediction 
of a quasi universal a\ parameter. However, a reduction of the experimental 
uncertainties to the percent level would be very desirable for obtaining a more 
conclusive picture. 

We start by considering ratios of non-leptonic decay rates that are related 
to each other by the replacement of a pseudoscalar meson by a vector meson. 
In the comparison of B — > Dir and B — > D*ir decays one is sensitive to the 



difference of the values of the two a\ parameters in (64) evaluated for a{ = 0. 
This difference is at most few times 10 -3 . Likewise, in the comparison of 
B — > Dir and B — > Dp decays one is sensitive to the difference in the light-cone 
distribution amplitudes of the pion and the p meson, which start at the second 
Gegenbauer moment a\. These effects are suppressed even more strongly. 
From the explicit expressions for the decay amplitudes in (|58|) it follows that 



T(B d ^D+7T-) _ (m 2 B - m 2 D ) 2 \q\ D „ { F (ml) 



T{B d ^D*+v-) 4m 2 B \q\ 3 D ^ \Ao(ml) 



ai(£>7r) 



ai(£»*?r) 



T(B d - 


-D+p-) 




f 2 


(F + (m 2 p ) 




ai {Dp) 


r{B d - 


-> D+TT-) 


(m% - m 2 D ) 2 \q\ D K 


f 2 

J TV 









(65) 

Using the experimental values for the branching ratios reported by the CLEO 
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CollaborationE] we find (taking into account a correlation between some sys- 
tematic errors in the second case) 



cn(Dir) 



ai(£»*7r) 
ai(Dp) 



A (ml) 



= 1.00 ±0.11. 



ai(L>7r) 



F+(m 2 ) 

; P X = 1.16±0.11 

F Q {m%) 



(66) 



Within errors, there is no evidence for any deviations from naive factorization. 

Our next-to-leading order results for the quantities a\{D^*'L) allow us 
to make theoretical predictions which are not restricted to ratios of hadronic 
decay rates. A particularly clean test of these predictions, which is essentially 
free of hadronic uncertainties, is obtained by relating the B d — > D^ + L^ decay 
rates to the differential semi-leptonic B d — > £>(*)+ l~v decay rate evaluated at 
q 2 = m 2 L . In this way the parameters |ai| can be measured directly 113. One 
obtains 



dT{B d -> D^+l~v)/dq 2 



= 6n 2 \V ud \ 2 f 2 L \ ai (D^L)\ 2 X 



-(*) 



(67) 

where X p = X* = 1 for a vector meson (because the production of the lepton 
pair via a V — A current in semi-leptonic decays is kinematically equivalent to 
that of a vector meson with momentum q) , whereas X n and X* deviate from 1 
only by {calculable) terms of order m 2 / m 2 B , which numerically are below the 
1% levellia. We emphasize that with our results for a\ given in (43) the above 
relation becomes a prediction based on first principles of QCD. This is to be 
contrasted with the usual interpretation of this formula, where a\ plays the 
role of a phenomenological parameter that is fitted from data. 

The most accurate tests of factorization employ the class-I processes B d — * 
D* + L~, because the differential semi-leptonic decay rate in B — > D* transi- 
tions has been measured as a function of q 2 with good accuracy. The results of 
such an analysis, performed using CLEO data, have been reported inE3. One 
finds 



R* = (1.13 ± 0.15) GeV 2 
R* = (2.94 ± 0.54) GeV 2 
(3.45 ± 0.69) GeV 2 



R ai 



|oi(D*tt)| = 1.08 ±0.07, 
\ ai (D*p)\ = 1.09 ±0.10, 
|ai(D*oi)| = 1.08 ±0.11. (68) 



This is consistent with our theoretical result in (j43|). In particular, the data 
show no evidence for large power corrections to our predictions obtained at 
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Table 3: Model-independent predictions for the branching ratios (in units of 10 -3 ) of class- 
I, non-leptonic —* L~ decays in the heavy-quark limit. All predictions are in units 

of (|oi |/1.05) 2 p-iThe last two columns show the experimental results reported by the CLEO 
Collaboration^, and by the Particle Data GroupEj. 



Decay mode 


Theory (HQL) 


CLEO data PDG98 


B d -> D+tt- 
B d -> D+R- 
B d -> D+p- 
B d -> D+R*- 
B d -> D+a^ 


3.27 
0.25 
7.64 
0.39 
7.76 
x[F + (0)/0.6] 2 


2.50 ±0.40 3.0 ±0.4 
7.89 ± 1.39 7.9 ±1.4 
8.34 ± 1.66 6.0 ±3.3 


B d -► D*+tt- 
B d -> D*+R- 
B d -> Z>*+p- 
£ d -> D*+R*- 
B d -> £*+ar 


3.05 
0.22 
7.59 
0.40 
8.53 
x[A (0)/0.6] 2 


2.34 ±0.32 2.8 ±0.2 
7.34 ± 1.00 6.7 ±3.3 
11.57 ±2.02 13.0 ±2.7 



leading order in Aqcd/^6- However, a further improvement in the experi- 
mental accuracy would be desirable in order to become sensitive to process- 
dependent, non-factorizable effects. 

7.4 Predictions for class-I decay amplitudes 

We now consider a larger set of class-I decays of the form B d — > D^ + L~ , all 
of which are governed by the transition operator (^) . In Table || we compare 
the QCD factorization predictions with experimental data. As previously we 
work in the heavy-quark limit, i.e. our predictions are model independent up 
to corrections suppressed by at least one power of Aqcd/ti&- The results show 
good agreement with experiment within errors, which are still rather large. 
(Note that we have not attempted to adjust the semi-leptonic form factors 
F + (0) and Aq(0) so as to obtain a best fit to the data.) 

We take the observation that the experimental data on class-I decays into 
heavy-light final states show good agreement with our predictions obtained in 
the heavy-quark limit as evidence that in these decays there are no unexpect- 
edly large power corrections. We will now address the important question of the 
size of power corrections theoretically. To this end we provide rough estimates 
of two sources of power-suppressed effects: weak annihilation and spectator in- 
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teractions. We stress that, at present, a complete account of power corrections 
to the heavy-quark limit cannot be performed in a systematic way, since these 
effects are not dominated by hard gluon exchange. In other words, factoriza- 
tion breaks down beyond leading power, and there are other sources of power 
corrections, such as contributions from higher Fock states, which we will not 
address here. We believe that the estimates presented below are nevertheless 
instructive. 

To obtain an estimate of power corrections we adopt the following, heuris- 
tic procedure. We treat the charm quark as light compared to the large scale 
provided by the mass of the decaying b quark (m c -C and m c fixed as 
mi, — ► oo) and use a light-cone projection similar to that of the pion also 
for the D meson. In addition, we assume that m c is still large compared to 
Aqcd- We implement this by using a highly asymmetric D-meson wave func- 
tion, which is strongly peaked at a light-quark momentum fraction of order 
Aqcd/w_d. This guarantees correct power counting for the heavy-light final 



states we are interested in. As discussed in Sect. 5.2, there are four annihila- 
tion diagrams with a single gluon exchange (see Fig. || (a)-(d)). The first two 
diagrams are "factorizable" and their contributions vanish because of current 
conservation in the limit m c — > 0. For non-zero m c they therefore carry an 
additional suppression factor m 2 D /m 2 B « 0.1. Moreover, their contributions to 
the decay amplitude are suppressed by small Wilson coefficients. Diagrams (a) 
and (b) can therefore safely be neglected. From the non-factorizable diagrams 
(c) and (d) in Fig. || the one with the gluon attached to the b quark turns out 
to be strongly suppressed numerically, giving a contribution of less than 1% 
of the leading class-I amplitude. We are thus left with diagram (d), in which 
the gluon couples to the light quark in the B meson. This mechanism gives 
the dominant annihilation contribution. (Note that by deforming the light 
spectator-quark line one can redraw this diagram in such a way that it can be 
interpreted as a final-state rescattering process.) 

Adopting a common notation, we parameterize the annihilation contribu- 
tion to the Bd — > D + it~ decay amplitude in terms of a (power-suppressed) am- 
plitude A such that A(Bd — ► D + ir~) = T + A, where T is the "tree topology", 
which contains the dominant factorizable contribution. A straightforward cal- 
culation using the approximations discussed above shows that the contribution 
of diagram (d) is (to leading order) independent of the momentum fraction £ 
of the light quark inside the B meson: 

A ~ UJdJb / du / dv ~ 3f w f D jB / dv . (69) 

J u J v z J v z 

The B- meson wave function simply integrates to fs, and the integral over the 
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pion distribution amplitude can be performed using the asymptotic form of the 



wave function. We take $d(v) in the form of (62) with the coefficients af = 0.8 
and a® — 0.4 (af = 0, i > 2). With this ansatz is strongly peaked at 

v ~ Aqcd/W-D. The integral over $d(v) in ( |69| ) is divergent at v = 1, and 
we regulate it by introducing a cut-off such that v < 1 — A/m B with A « 
0.3 GeV. Then J" $d (u)/v 2 w 34. Evidently, the proper value of A is largely 
unknown, and our estimate will be correspondingly uncertain. Nevertheless, 
this exercise will give us an idea of the magnitude of the effect. For the ratio of 
the annihilation amplitude to the leading, factorizable contribution we obtain 

A^a s C ++ C. f DfB f dv ijM* 0M . (to) 



T 3 2C+ + C- F o (0)m 



We have evaluated the Wilson coefficients at /i = nib and used fn = 0.2 GeV, 
f B = 0.18 GeV, F (0) = 0.6, and a s — 0.4. This value of the strong coupling- 
constant reflects that the typical virtuality of the gluon propagator in the 
annihilation graph is of order Aqcd^b- We conclude that the annihilation 
contribution is a correction of a few percent, which is what one would expect 
for a generic power correction to the heavy-quark limit. Taking into account 
that f B ~ A Q cD(AQCD/m i3 ) 1 / 2 , F (0) ~ (A QC D/m B ) 3 / 2 and f D ~ A QCD , we 
observe that in the heavy-quark limit the ratio AjT indeed scales as Aqcd /mb, 
exhibiting the expected linear power suppression. (Recall that we consider the 
D meson as a light meson for this heuristic analysis of power corrections.) 

Using the same approach, we may also derive a numerical estimate for 
the non-factorizable spectator interaction in Bd — > D + ir~ decays, discussed in 



Sect. 5.1. We find 



^spec 27TQ S C+ - C- f D f B m B <f> D {v) , , 

m — ~ -^r~ on _l n f 7n\ 2 x / dv — ^--0.03, (71 
Ti ead 3 26+ + 6_ F o (0)m B A 



B 



where the hadronic parameter \ B = O(Aqcd) is defined as ^(d^/O ^s(0 = 
m B /X B . For the numerical estimate we have assumed that X B 0.3 GeV. 
With the same model for <£>£>(«) as above we have J dv $>d(v)/v w 6.6, where 
the integral is now convergent. The result (|7l|) exhibits again the expected 
power suppression in the heavy-quark limit, and the numerical size of the 
effect is at the few percent level. 

We conclude from this discussion that the typical size of power corrections 
to the heavy-quark limit in class-I decays of B mesons into heavy-light final 
states is at the level of 10% or less, and thus our prediction for the near 
universality of the parameters a\ governing these decay modes appears robust. 
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7.5 Remarks on class-II and class-Ill decay amplitudes 

In the class-I decays Bd — * D^*' + L~ , the flavour quantum numbers of the 
final-state mesons ensure that only the light meson L can be produced by the 
(du) current contained in the operators of the effective weak Hamiltonian in 
(|l5|). The QCD factorization formula then predicts that the corresponding 
decay amplitudes are factorizable in the heavy-quark limit. The formula also 
predicts that other topologies, in which the heavy charm meson would be 
created by a (cu) current, are power suppressed. To study these topologies 
we now consider decays with a neutral charm meson in the final state. In the 
class-II decays Bd — * D^°L° the only possible topology is to have the charm 
meson as the emission particle, whereas for the class-Ill decays B~ — > D^ L~ 
both final-state mesons can be the emission particle. The factorization formula 
predicts that in the heavy-quark limit class-II decay amplitudes are power 
suppressed with respect to the corresponding class-I amplitudes, whereas class- 
Ill amplitudes should be equal to the corresponding class-I amplitudes up to 
power corrections. 

It is convenient to introduce two common parameterizations of the decay 
amplitudes, one in terms of isospin amplitudes Ai/ 2 and A 3 / 2 referring to 
the isospin of the final-state particles, and one in terms of flavour topologies 
(T for "tree topology", C for "colour suppressed tree topology", and A for 
"annihilation topology" ) . Taking the decays B — > Dtt as an example, we have 

A(B d -> D+7T- ) = J^A 3/2 + \/|^i/2 - T + A , 
V2A(B d - D°n°) = \[\a z/2 - y^A 1/2 = C - A, 

A(B~ -> D a ir- ) = V3A 3/2 = T + C . (72) 

A similar decomposition holds for the other B — > D^L decay modes. Isospin 
symmetry of the strong interactions implies that the class-Ill amplitude is 
a linear combination of the class-I and class-II amplitudes. In other words, 
there are only two independent amplitudes, which can be taken to be Ai/ 2 
and A 3 / 2 , or (T + A) and (C — A). These amplitudes are complex due to 
strong-interaction phases from final-state interactions. Only the relative phase 
of the two independent amplitudes is an observable. We define 6 to be the 
relative phase of A 1 / 2 and A z / 2 , and Stc the relative phase of (T + A) and 
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Table 4: CLEO data BO on the branching ratios for the decays B — > D^L in units of 
10 -3 . Upper limits are at 90% confidence level. See text for the definition of the quantities 
<5 and TZ. 





B^Dix B^Dp B^D*n B^D*p 


Class-I (£>«+£") 
Class-II (L>M L ) 
Class-Ill (D(*y°L-) 


2.50 ± 0.40 7.89 ± 1.39 2.34 ±0.32 7.34 ±1.00 

< 0.12 < 0.39 < 0.44 < 0.56 
4.73 ± 0.44 9.20 ±1.11 3.92 ± 0.63 12.77±1.94 


5 

n 


< 22° < 30° < 57° < 31° 
1.34 ±0.13 1.05 ±0.12 1.26 ±0.14 1.28 ±0.13 



(C — A). The QCD factorization formula implies that 

A/2 



V2A 3/2 

C-A 
T + A 



= 1 + 0(A QC D/m 6 ) , 5 = 0(A QC d/toe,) . 



(9(A QCD /m fc ) , S TC = 0(1) . (73) 



In the remainder of this section, we will explore to what extent these predictions 
are supported by data. 

In Table || we show the experimental results for the_yarious B — > D^L 
branching ratios reported by the CLEO Collaboration EM3. We first note that 
no evidence has been seen for any of the class-II decays, in accordance with 
our prediction that these decays are suppressed with respect to the class-I 
modes. Below we will investigate in more detail how this suppression is real- 
ized. The fourth line in the table shows upper limits on the strong-interaction 
phase difference 8 between the two isospin amplitudes, which follow from the 



relation 



sin 2 * < I T -^l Bl '^ - D% °^ (74) 

The strongest bound arises in the decays B — > Dn, where the strong-interaction 
phase is bound to be less than 22°. This confirms our prediction that the phase 
6 is suppressed in the heavy-quark limit. 

Let us now study the suppression of the class-II amplitudes in more detail. 



We have already mentioned in Sect. 6.4 that the observed smallness of these 



amplitudes is more a reflection of colour suppression than power suppression. 
This is already apparent in the naive factorization approximation, because the 
appropriate ratios of meson decay constants and semi-leptonic form factors 
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exhibit large deviations from their expected scaling laws in the heavy-quark 
limit, as shown in (|57|). Indeed, it is obvious from Table [| that there are 
significant differences between the class-I and class-Ill amplitudes, indicating 
that some power-suppressed contributions are not negligible. In the last row 
of the table we show the experimental values of the quantity 



K 



A{B- 



D(*)°L~ 



A(B d -> Di*)+L- 



I T(B d ) Br(B- -> D(*)°L-) 
t{B~) Bi(B d -> D(*)+L-) 



(75) 



which parameterizes the magnitude of power-suppressed effects at the level 
of the decay amplitudes. If we ignore the decays B — > D* p with two vector 
mesons in the final state, which are more complicated because of the presence 
of different helicity amplitudes, then the ratio 1Z is given by 



K = 



1 



C-A 



T + A 



0.2 

l + x — 

Ol 



(76) 



where a\ are the QCD parameters entering the transition operator in ( [40| ) , and 
a-2 = ^l^j ~ C+ — " non_ f ac t° r i z abl e corrections" (77) 

jLi 1 V f £t J V f 



are the corresponding parameters describing the deviations from naive factor- 
ization in the class-II decays E3. All quantities in (|76| ) depend on the nature of 
the final-state mesons. In particular, the parameters 



x(Dtt) 



x{Dp) 



x{D*ir) 



(m? B -ml)f D F B ^(m 2 D ) 
(m 

f D A a 



'I 



0.9, 



f p F^ D (mi p ) 

f D ,Ff^(m 2 D ,) 
UA^ D '{ml) 



0.5, 



0.9, 



(78) 



account for the ratios of decay constants and form factors entering in the naive 
factorization approximation. For the numerical estimates we have assumed 
that the ratios of heavy-to-light over heavy-to-heavy form factors are approxi- 
mately equal to 0.5, and we have taken jo = 0.2 GeV and fr>* — 0.23 GeV for 
the charm meson decay constants. Note that in ([76]) it is the quantities x that 
are formally power suppressed ~ (Aqcd/w{,) 2 in the heavy-quark limit, not 
the ratios a^jax. For the final states containing a pion the power suppression 
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is clearly not operative, mainly due to the fact that the pion decay constant f n 
is much smaller than the quantity (/d-\/«i_d) 2 / 3 ~ 0.42 GeV. To reproduce the 
experimental values of the ratios 1Z shown in Table ^ requires values of a%/a\ 
of order 0.1-0.4 (with large uncertainties), which is consistent with the fact 
that these ratios are of order 1/N C in the large- N c limit, i.e. they are colour 
suppressed. 

The QCD factorization formula (|3[) allows us to compute the coefficients ai 
in the heavy-quark limit, but it does not allow us to compute the corresponding 
parameters 02 in class-II decays. The reason is that in class-II decays the 
emission particle is a heavy charm meson, and hence the mechanism of colour 
transparency, which was essential for the proof of factorization, is not operative. 
For a rough estimate of ai in B — > ttD decays we consider as previously the 
limit in which the charm meson is treated as a light meson (m c <C nib), however 
with a highly asymmetric distribution amplitude. In this limit wej-jcan adapt 
our results for the class-II amplitude in B — > tttt decays derived inH, with the 
only modification that the hard-scattering kernels must be generalized to the 
case where the leading-twist light-cone distribution amplitude of the emission 
meson is not symmetric. We find that 

2N C +VW 2N C VW 



( ' F [d+(M) + C-Qi)] (~61n£j + /, + /„ ) . (7fh 



4tt 2N c 



where 



// = / dv 
Jo 



\n 2 v + In v H 6 + m{2 In v - 3) + 0(v) 

3 



12tt 2 ,Uf B m B f , $d{v) 



fl1 = IV" FB^ (m 2 )m 2 ~~\ / dv — — ■ ( 80 ) 

l\ c t Q (m D )m B ab J v 

The contribution from /// describes the hard, non-factorizable spectator in- 
teraction. Note that this term involves J dv$>ij(v)/v, which can be sizeable 
but remains constant in the heavy-quark limit implied here (mj — > 00 with 
m c constant). Using the same numerical inputs as previously, we find that 
fu ~ 13 and fi ~ —1 — I9i. In writing the hard-scattering kernel for fj 
we have only kept the leading terms in v, which is justified because of the 
strongly asymmetric shape of <&£>(«). Note the large imaginary part arising 
from the "non-factorizable" vertex corrections with a gluon exchange between 
the final-state quarks. Combining all contributions, and taking u, — nib for the 
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rcnormalization scale, we find 



a 2 « 0.25e~ l41 ° , (81) 

which is significantly larger in magnitude than the leading-order result aJj ss 
0.12 corresponding to naive factorization. We hasten to add that our estimate 
( |Sl| ) should not be taken too seriously, since it is most likely not a good ap- 
proximation to treat the charm meson as a light meson. Nevertheless, it is 
remarkable that in this idealized limit one obtains indeed a very significant 
correction to naive factorization, which gives the right order of magnitude for 
the modulus of a 2 and, at the same time, a large strong-interaction phase. 



For completeness, we note that the value for a 2 in (81) would imply a strong- 
interaction phase difference S « 10° between the two isospin amplitudes A 1 / 2 
and A 3 / 2 in B — > Dn decays, and hence is not in conflict with the experimental 
upper bound on the phase S given in Table ||. The phase 5tc i 011 the other 
hand, is to leading order simply given by the phase of a 2 and is indeed large, 
in accordance with (|7E 



8 Conclusion 

With the recent commissioning of the B factories and the planned emphasis 
on heavy-flavour physics in future collider experiments, the role of B decays in 
providing fundamental tests of the Standard Model and potential signatures 
of new physics will continue to grow. In many cases the principal source of 
systematic uncertainty is a theoretical one, namely our inability to quantify 
the non-perturbative QCD effects present in these decays. This is true, in 
particular, for almost all measurements of CP violation at the B factories. 

In these lectures, I have reviewed a rigorous framework for the evaluation 
of strong- interaction effects for a large class of exclusive, two-body non-leptonic 
decays of B mesons. The main result is contained in the factorization formula 
, which expresses the amplitudes for these decays in terms of experimentally 
measurable semi-leptonic form factors, light-cone distribution amplitudes, and 
hard-scattering functions that are calculable in perturbative QCD. For the first 
time, therefore, we have a well founded field-theoretic basis for phenomenolog- 
ical studies of exclusive hadronic B decays, and a formal justification for the 
ideas of factorization. For simplicity, I have mainly focused on B — > Dtt decays 
here. A detailed discussion of B decays into two light mesons will be presented 
in a forthcoming paper 113. 

We hope that the factorization formula (J3J) will form the basis for future 
studies of non-leptonic two-body decays of B mesons. Before, however, a fair 
amount of conceptual work remains to be completed. In particular, it will be 
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important to investigate better the limitations on the numerical precision of the 
factorization formula, which is valid in the formal heavy-quark limit. We have 
discussed some preliminary estimates of power-suppressed effects in the present 
work, but a more complete analysis would be desirable. In particular, for rare 
B decays into two light mesons it will be important to understand the rolp pf 
chirally-enhanced power corrections and weak annihilation contributions 
For these decays, there are also still large uncertainties associated with the 
description of the hard spectator interactions. 

Theoretical investigations along these lines should be pursued with vigor. 
We are confident that, ultimately, this research will result in a theory of non- 
lcptonic B decays, which should be as useful for this area of heavy-flavour 
physics as the large-Wf, limit and heavy-quark effective theory were for the 
phenomenology of semi-lcptonic decays. 
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